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1. EXECUTIVE SUMMARY 

A team consisting of collaborators from Tennessee Technological University, MSFC, BD Systems and 
the University of Delaware (herein called the TTU team) conducted specific research and development 
activities in MXER tether systems during the base period of 5-15-04 through 3-24-06 under contract 
(numbers NNM04AB 13C). The team addressed two primary topics related to the MXER tether 
system; development of validated, high-fidelity dynamic models of an elastic rotating tether and 
development of feasible mechanisms to enable reliable rendezvous and capture. The following report 
will describe in detail the activities that were performed during the base period of this cycle 2 MXER 
tether activity, and will summarize the results of this funded activity. A brief summary of the 
significant outcomes of this phase I activity are provided here. 

MXER Tether Capture Mechanisms: 

• Delivery of Quad-Trap , a robust capture mechanism proposed, developed, tested and 
demonstrated with high degree of feasibility 

• High-fidelity dynamic model of the capture mechanism created and validated through 
comparisons with hardware 

• Autonomous capture demonstrated using the Quad-trap mechanism 

• Significant number (40+) of unique capture concepts proposed and evaluated 

• Trade study tool established for evaluating proposed capture concepts 

• Hierarchal (based on technical aspects) categorization tool developing for organizing capture 
concepts 

• A series of other feasible capture mechanisms developed and described 

• Passive operation (Gravity-gradient actuation) demonstrated in capture mechanism 

• Test facility developed to enable kinematically & dynamically-similar testing of capture 
systems 

MXER Tether Dynamic Modeling 

• Detailed development of high-fidelity elastic tether dynamic models provided through multiple 
formulations 

• Detailed development of specific algorithms for efficient implementation of tether dynamic 
models 

• Analysis of error window at capture performed based on Monte Carlo simulation with tether 
dynamic model used to guide subsequent capture mechanism synthesis 

• Validation procedure implemented and conducted based on cross-model of verification 

• Observed several, potentially significant, non-linear effects in tether dynamics not currently 
reported in the literature or in assumed-modes methods: 

o Coupling between axial and transverse modes resulting in resonance or beating 
o System response dependent on initial conditions in a non-linear fashion 
o Resonant conditions possible in all tether configurations 

• Developed a series of system response results that can be used to validate other tether 
algorithms. 

• One (1) paper ready to submit to ISPT export control in July, 2004 for publication, additional 
papers to be submitted in August 2005 to ISPT export control for publication. 
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MXER Tether System Components 

• Spatial 3 dof mechanism proposed, developed and demonstrated mechanism to achieve 
efficient sun-tracking by tether-node-based solar arrays 

• Three prototype models with real-time tracking developed and delivered 

• Dynamic models and detailed design presented for practical implementation 

• One prototype model developed and tested for possible attitude control system 

2. OVERVIEW OF MXER SYSTEM 

Momentum-exchange/electrodynamic reboost (MXER) tether systems may enable high-energy 
missions to the Moon, Mars and beyond by serving as an “upper stage in space”. Existing rockets that 
use a MXER tether station could double their capability to launch communications satellites and help 
improve US competitiveness. A MXER tether station would boost spacecraft from low Earth orbit to a 
high-energy orbit quickly, like a high-thrust rocket. Then, using the same principles that make an 
electric motor work, it would slowly rebuild its orbital momentum by pushing against the Earth's 
magnetic field— without using any propellant. While many significant challenges exist in designing a 
feasible MXER tether system, the proposal of a passive capture mechanism [3] provided the impetus 
necessary to consider MXER tethers for future in-space transportation. With this approach, the 
significant challenges in developing a momentum-exchange/electrodynamic reboost tether system are 
in the ability to predict the behavior of the system in a manner to control tether flight and coordinate a 
successful rendezvous between tether and payload. Accurate tether propagation, combined with a 
successful capture mechanism will provide two of the keys in making the technology feasible for 
future application. 


3. PROJECT STATUS 

According to task statement and deliverables 
TTU performing organization 


Task 

# 

Description 

Status 

Date 

Completed 

Source Materials 

Tether Dynamic Modeling and Algorithm Activity 



1.1 

Develop the 
MXER tether 
system dynamic 
models and 
algorithms to be 
used within the 
complete tether 
system model 

Four analytical models developed to 
consider tether dynamics 

1 . Continuum model using inertial- 
frame-based coordinates and finite 
difference (FD) solver 

2. Continuum model using relative- 
orbit-based coordinates and finite 
difference (FD) solver 

3. Continuum model applicable to 
multi-noded tether system, inertial- 
based coordinates, finite difference 
solver 

Models to include axial and bending 
stiffness, various damping types. 

Three classes of computational ODE 
algorithms are implemented: 

1 . Runge Kutta (single-point explicit) 

12-1-05 

Model development 
described in section 
B 1, results of these 
models can be seen in 
section B 2, B 3 and 
Appendix B. 
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Task 

# 


Description 


Status 


Date 

Completed 


Source Materials 




2. Adams-Bashforth Moulton (multi- 
point explicit) 

3. Gear Methods (multi-Point implicit) 
For the implicit techniques, the ODE 
system gradient is evaluated for use in a 
Newton-type search 



1.2 

Validation of 
tether system 
models, with 
extended 
implementation 

Multi-stage process of model validation 
defined through the “Model Validation 
Matrix”. This includes: 

1) compare with closed-form solution 
(Wave equation) 

2) Spatial/Temporal convergence trials 

3) Fimiting MXER configs with CFS 

4) Compare with in-Fab experiments 

5) Comparison of reference MXER 
applications 

This validation will spans TTU team 
model development activities and is 
available for other participants. 

2-15-06 

(Validation 

material 

contributed 

by all team 

members) 

Validation activity 
described in detail in 
section B 2 with 
results included. 

1.3 

Basic 

implementation 

and 

demonstration in 
selected software 
platform. 

Analytical Models 1 and 4 implemented 
and demonstrated in Matlab using three 
classes of ODE solvers. 

11-15-05 

Implementation given 
in section B 3 

1.4 

Basic application 
of the dynamic 
model as a 
simulation, 
analysis and 
design tool. 

Initial used of tether model as a 
simulation, analysis and design tool 
performed in evaluation of tether 
capture window, evaluated as a 
stochastic predictor process using 
Monte-Carlo simulation. 

12-15-05 

Demonstration as use 
as a design tool given 
in section A 1 and 
Appendix B. 

1.5 

Develop initial 
condition and 
boundary 
condition models 

Models for initial conditions developed, 
providing closed-form solution for 
initial rotational state. Results delivered 
in algorithmic format along with PDE 
and its solution. 

2-25-05 

Material presented in 
Section B 5 

1.6 

Participate in 
integration of 
tether dynamic 
model and 
development of 
the MXER tether 
system design 
and analysis tool. 

Integration activities ongoing 
throughout Base period extension. 

3-15-06 


1.7 

Final reporting 
and presentation 
of year one 
results shall be 
prepared in 
accordance with 

Mid-term report delivered 10-13-04 
Base-year Final Report deliver 7-15-05 
Project final report, delivery date of 5- 
15-06 (draft 3-31-06) 

5-15-06 
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Task 

# 


Description 


Status 


Date 

Completed 


Source Materials 



DRD 1040MA- 
001. 




Capture Mechanism Development Activity: 

2.1 

Enumerate 
objective 
functions, figures 
of merit and 
associated tether 
behavior that will 
drive the design 
and development 
of capture system 
concepts. 

Task statement and design objectives 
prepared to guide capture system design 
and development process. These are 
demonstrated in part as evaluation 
points in the capture-concept trade 
study. One significant capture objective 
was definition of a probabilistic capture 
window based on tether flight 
dynamics, as well as clear definition of 
the kinematic and dynamic nature of the 
tether capture process. 

10-15-05 

Section A 1 

2.2 

Define an array of 
capture system 
concepts. 

Initial set of capture concepts developed 
and proposed. The final set includes 40 
+ concepts. Each concept or class of 
concepts in the set is defined using the 
quad-chart format. Organization of all 
capture concepts is provided through 
the development of the capture concept 
genealogy and naming system. The 
genealogy and the capture concept 
description library are hyperlinked for 
easy review. 

6-15-05 

Section A 2 

2.3 

Conduct a capture 
system trade 
study. 

A complete trade study on the capture 
concepts proposed was performed. 
From this, four concepts were selected 
for additional analysis. 

1-15-05 

Section A 3 

2.4 

Select and 
develop basic 
designs for most 
promising capture 
concepts selected 
from Task 2.3. 

Four basic designs were selected from 
the first stage capture concept trade 
study; 

1. Modified IRIS (Quad-Trap) 
(CCM.T074b.P03x. A) 

2. A modified Turkey Trap 
(CCM.T03.P05b.A) 

3. Modified V-Gripper 
(CCM.T03d.P04.P) 

4. Tethered Goal-post flier 
(CCM.T03.P02b.P) 

An engineering analysis and simple 
design of each concept was performed. 
These results were used to down select 
to the top capture concept candidate. 
Advanced models and further trade 
study analysis performed based on this 
work. 

2-15-05 

Section A 4.1 

through section A 4.6 

2.4.1 

Evaluation of 3 
dof spatial joint 
mechanism for 

1. Two demonstration prototype 
mechanisms fabricated and delivered, 
one capable of real-time control by 

2-15-06 

Section A 4.7 

Simulations 

delivered: 


4 



Task Description 
# 


Status 


Date 

Completed 


Source Materials 



MXER and other 
in-space 
propulsion 
concepts 

operator through joystick input. 

2. Kinematic models of mechanism 
constructed and made available 

3. Issues in design for as an in-space 
mechanism considered 

4. Alternative applications related to in- 
space activities considered and 
addressed: Kinematic analysis and 
simulation of various in-space uses with 
supporting data provided (mod. 4). 
Develop and demonstrated prototype 
system based on new design with 
internal, pressurized fluid passageways, 
tested as a mounting system for attitude 
control thrusters 


2 simulations for 
Bigelow architecture, 

2 general mechanism 
simulations (thruster 
applications), n 
simulations for 
multiple arrays. 

3 prototypes 
delivered, one 
integrated with 
thruster for testing. 

2.5 

Perform testing 
and demonstration 
on proof-of- 
concept 
prototypes of 
capture system 
concepts 

developed in task 
2.4. 

Testing and demonstration initiated 

through the following activities: 

1 . Design, construction and testing of 
the payload launch system (PLS) 
and test facilities - 11-15-04 

2. Design and construction of the 
payload - 12-13-04 

3. Design, construction, testing and 
demonstration of the capture 
mechanism. 4-15-05 

4. Testing of LM mechanisms (8-05) 

5. Advanced dynamic testing of the 
Quad-Trap (2-06) 

2-15-06 

Section A 5 

2.5.1 

Redesign of 
capture concepts 
from year one 
with advanced 
testing and 
demonstration of 
capture system 
prototype 

Advanced developments of the capture 
mechanism were addressed in the 
following activities: 

1 . Detailed dynamic testing and analysis 
of the capture mechanism, leading to 
advanced model and system 
understanding (2-15-06) 

2. Continued work in autonomous 
operation of the capture mechanism 
with two alternative autonomous 
capture systems demonstrated, 3-24-06 

3. Testing of alternative capture 
mechanisms 8-15-05 

3-24-06 

Section A 5.8 

2.5.2 

Investigate issues 
resulting from 
capture concepts 
developments in 
year one 

The issue associated with slender 
booms was identified as one of the 
primary areas for advanced 
investigation on year-one capture 
mechanism development. Slender 
booms were addressed through 
activities in modeling, analysis, 
characterizing behavior and issues of 
practical implementation as associated 
with the Quad-trap mechanism. 

3-29-06 

Section A 4.1, 4.2 
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Task 

Description 

Status 

Date 

Source Materials 

# 



Completed 


2.6 

Final reporting 
and presentation 
of year one results 
shall be prepared 
in accordance 
with DRD 
1040MA-001. 

Mid-term report delivered 10-13-04 
Final report delivered 11-9-06. 




4. TTU TEAM AND PRIMARY PARTNER CONTRIBUTIONS 


Name 

PI/Co-I 

Org 


Work Commitments 

Stephen Canfield 

PI 

TTU 

• 

• 

• 

• 

Overall project manager, organize activities to 
meet technical objectives and deliverables 
Lead effort in capture system trade study 
Participate in developing specific capture 
concepts 

Participate in capture model development 

Jamie Beard 


TTU 

• 

John Peddieson 

Co-I 

TTU 

• Participate in tether dynamic model 
development and validation 

Joe Richardson 

Co-I 

TTU 

• Participate in capture system trade study 

• Participate in tether model development 

Jamie Beard 

Research 

TTU 

• Computational implementation, capture model 
development and testing 

Mike Renfro 

Research 

TTU 

• httD://mxer.tntech.edu maintenance 

Seth Knight 

G.S. 

TTU 

• Capture system development, prototype testing 

Marshall Norris 

G.S. 

TTU 

• Capture system development, quad trap model 
(ADAMS) 

Daniel Chlarson 

G.S. 

TTU 

• Prep Facilities, testing, report editing 

Neil White 

G.S. 

(hourly) 

TTU 

• Developed PLS, fabrication 

Sunil Agrawal 

Co-I 

UD 

• Participate in developing specific capture 
concepts 

• Participate in dynamic model validation work 

John Glaese 

Co-I 

BDSys 

• 

• 

Participate in developing specific capture 
concepts 

Participate in capture model development 
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5. REPORT SUMMARY 


The remainder of this report describes in detail the activities that were performed during the base year 
to develop capture mechanisms and dynamic models for the MXER tether system. The report is 
divided into two sections; Section A covers development of the capture mechanism while Section B 
covers the dynamic model development and validation activities. Each section proceeds in a 
sequential order according to the tasks statement, i.e., section A 1.0, Derive Capture System Objectives 
describes the work conducted to complete task statement 1 on the capture mechanism, and so-forth. 
Two additional sections are included at the end of the report providing an acronym list, a derivation of 
the Error Space prediction for capture, and CAD model drawing of the Quad Trap Mechanism 
prototype. A results and conclusions section is found near the end of the report. 
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1.0 DERIVE OBJECTIVES 

The first task derives the objectives that will drive the design and development of capture system 
concepts. 

1.1 Evaluation Criteria for the Capture Mechanism 

The capture mechanism design and development process is guided by a task statement and design 
objectives. These were prepared from a combination of system requirements, kinematic and dynamic 
description of the capture event, and an evaluation of the expected error window in capture. The 
specific outcomes of these considerations are demonstrated in part as evaluation points in the capture- 
concept trade study. The evaluation criteria defined for conducting the trade study are given in Table 
1 . 


Table 1: Capture Mechanism Evaluation Criteria 


Criterion 

Definition 

Low Mass Tether 

Physical mass of the portion of the capture device which rides on the 
tether alone 

Low Mass Payload 

Physical Mass of the portion of the capture device which rides on the 
payload 

Spatial Expanse 

Volume over which the capture mechanism can effect a successful 
capture. 

DoF 

Out of plane 

Simplicity 

Overall simplicity of the design including number of parts, complexity 
of mechanism, active/ passive and complexity of programming 

Attempts Per 
Payload 

The number of attempts at a successful capture expected from a single 
payload launch 

Number of Uses 

The number of times a capture mechanism can be expected to boost a 
payload over the lifetime of the capture mechanism 

Current TRL 

Technical readiness of the proposed design at the current time 

Accuracy of 
Release 

Ability to effect the release of the payload at the time and in the 
direction desired 

Positive Grasp 

Ability of the mechanism to firmly latch the payload and tether 
masses. 

Dynamic Effects of 
Capture 

Any extraneous or unpredicted movement of the tether craft caused by 
the capture operation 

Dynamic Effects of 
Release 

Any extraneous or unpredicted movement of the tether craft caused by 
the release operation 

Reliability 

Probability of catch within the capture volume of the device 

Minimal Tether 
Considerations 

Any additional hardware that may need to be added to the tether craft 
in order to mate to the capture mechanism 

Minimal Payload 
Considerations 

Any additional hardware that may need to be added to the payload in 
order to mate to the capture mechanism 


1.2 Kinematic Description of the Capture Event 

The kinematic description of the capture event is described in the series of following figures, which are 
based on relative motions of the tether tip and payload prior to capture. The following three figures 
indicate the relative position of the tether with respect to the payload prior to capture ( x is local 
horizontal while y is local vertical) as well as the velocities near point of capture. Note that these are 
based on models that assume ideal, rigid-body motion of the tether-tip only. 
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Position at capture 



Figure 1.2-1: Relative Position between the Payload and Tether Tip Around Capture 


Position at capture 



Figure 1.2-2: Displacement Limits over 10 Seconds Surrounding Capture 
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Velocity at capture 



Figure 1.2-3: Relative Velocity between the Payload and Tether Tip Around Capture 

Finally, the tether tip motion relative to the payload (observer located on the payload) is described in a 
single graph (x vs. y) as shown in the figure below. This figure includes an estimate of the expected 
error in tether position, which will result from the analysis of the following section. 



Figure 1.2-4: Schematic of the Capture Event 

(Motion relative to payload) 

One significant capture objective was definition of a probabilistic capture window based on tether 
flight dynamics, as well as clear definition of the kinematic and dynamic nature of the tether capture 
process. The derivation of this error window is explained in Appendix A. Appendix A contains a 
formal paper prepared for submission to AIAA journal transactions titled “Error Space Prediction and 
Dynamic Response of Passive Capture Mechanisms in Tether Momentum Exchange”. 
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2.0 DERIVE AN ARRAY OF POTENTIAL CAPTURE SYSTEM CONCEPTS 


A large number of capture system concepts were developed by the extended TTU contract team to 
form a basis for the eventual capture mechanism. This section will present an overview of the capture 
concept mechanism development cycle as an iterative process through multiple generations. General 
areas of consideration are grouped and tenned concept families, which in turn generate new members 
of ideas. In addition, concepts perform cross over functions with other families to form new concepts. 
This process will be viewed in total through a genealogy chart. This chart is a living document and is 
updated to keep up with concept development. The genealogy is shown on page 17. A few notes on 
this figure are made; 

1) This genealogy attempts in include concepts that originate with TTU, LM and MSFC. A very 
general means of indicating the originator is given with a color scheme, black, grey and red for 
the concept specifier. Note that in many cases, the concepts originated from multiple points, at 
the moment this chart indicates such concepts also with black (TTU color) 

2) A naming scheme is generated for the concepts that consist of multiple parts. This naming 
scheme has been designed with the following requirements; 

a. Name indicates the general form and implementation of the concept 

b. Name indicates if the concept is passive or active 

c. Name indicates the general mechanism family 

d. Names allow room for expansion or insertion of new concepts. 

The name consists of the following components as shown in table 2: 


Table 2: Naming Scheme for the Capture Mechanism Genealogy 


CCM 

.Txxyz 

.Pxxyz 

.P/A 

Capture 

Mechanism 

Mechanism 

Passive or Active 

concept 

mechanism 

group on 
Tether 

group on 
payload 



xx- numeric digits indicating 
family 

y-alphanumeric digits indicating 
subgroup 

z -alphanumeric digits 
indication sub-subgroup 



The proposed capture system concepts were then categorized in the form of a genealogy chart (family 
tree) as shown on page 18. This method of categorization has the following attributes: 

• A number/naming scheme is developed that clearly identifies each concept and relates it to 
other concepts 

• All concepts are shown to be derived from a specific number of basic strategies of physical 
principles 

• The categorization scheme also differentiates components required on the tether and/or payload 

• Concepts developed by other groups (MSFC ISPT, LM) are included in this scheme 

• The genealogy chart is used to spawn new ideas. 
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3.0 CAPTURE SYSTEM TRADE STUDY 

This section will describe the approach used to perform a trade study on all capture system 
concepts. The trade study was used to identify the capture concepts that would most likely be 
successful at an early stage in the design process. In this way the trade study would prevent 
wasting limited resources on capture concepts which were not likely to succeed. The AHP 
(Analytical Hierarchy Process) was chosen to conduct this study in part because of its history in 
application to similar type engineering problems and because it can be applied to a large number 
of concepts quickly while these concepts are still at a very early stage of design. AHP also 
provided a graduated score for the trade study which would identify second and third place 
candidates to be tested. 

3.1 Brief Description of AHP 

AHP or Analytical Hierarchy Process is a means of objectively determining which option from a 
group of options is the most preferred for a particular set of criteria. This process is frequently 
used in business environments to make decisions where inputs are unclear or abstract. It is this 
ability to make an abstract decision objective which makes the AHP tool ideal for selecting a 
capture concept to build and test early in the design process. 

The tool is used as follows. First a set of design criteria are agreed upon. These criteria are 
specifications or a “wish list” for the resulting design. Each criterion is then pair wise compared 
to all other criteria to detennine a weighting for each. A set of choices or designs is then 
identified. Each choice is then pair wise compared to all other choices on the basis of the 
aforementioned criteria to weigh each choice on the basis of each individual criteria. A 
mathematical algorithm is then applied which outputs a final score for each design telling the 
user which design is most preferred. 

3.2 Overview of Trade Study Criteria 

Table 3 below provides and over and brief explanation of the trade study criterion chosen. 
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Table 3: Trade Study Criteria 


Criterion 

Low Mass Tether 

Low Mass Payload 
Spatial Expanse 
Simplicity 

Attempts Per Payload 
Number of Uses 

Current TRL 

Release mechanism 
Positive Grasp 
Dynamic Effects of Capture 

Dynamic Effects of Release 

Reliablity 

Minimal Tether Considerations 


Minimal Payload Considerations 


cost 

Economic impact 
Ground operations 

Sensors 


Definition 

Physical mass of the portion of the 
capture device which rides on the tether 
alone 

Physical Mass of the portion of the 
capture device which rides on the 
payload 

Volume over which the capture 
mechanism can effect a successful 
capture. 

Overall simplicity of the design including 
number of parts, complexity of 
mechanism, active/ passive and 
complexity of programming 
The number of attempts at a successful 
capture expected from a single payload 
launch 

The number of times a capture 
mechanism can be expected to boost a 
payload over the lifetime of the capture 
mechanism 

Technical readiness of the proposed 
design at the current time 
Ability of mechanism to cleanly release 
the payload. 

Ability of the mechanism to firmly latch 
the payload and tether masses. 

Any exteraneous or unpredicted 
movement of the tether craft caused by 
the capture operation 
Any exteraneous or unpredicted 
movement of the tether craft caused by 
the release operation 
Probability of device working as desired 
(measured as number of independent 
steps that need to occur) 

Any additional hardware that may need 
to be added to the tether craft in order to 
mate to the capture mechanism 

Any additional hardware that may need 
to be added to the payload in order to 
mate to the capture mechanism 

Economic impact for satellite provider 
Level of control/monitor needed on 
ground to make capture work 
Number and type of sensors needed 
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3.3 Results from the Trade Study 

The inputs into the AHP algorithm and the resulting scores are shown on Tables 4 and 5 below. 
However, for the convenience of the reader the following graph and tables are provided to 
illustrate the results of the AHP analysis. 


AHP Final Score* 

(Invert Values, Smaller is Better) 



□ Series 1 


Figure 3.3-1: Summary of the AHP Analysis and Final Scores 
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Mass 

Current 

TRL 

Probability 
of Capture 

Reliability 

Dynamic 
effects of 
capture and 
release 

Minimal 

tether/paylo 

ad 

considerati 

ons 

Const. Cost 

Mass 

1 

0.3333333 

0.1111111 

0.2 

0.3333333 

0.3333333 

7 

Current TRL 

3 

1 

0.1428571 

0.2 

0.3333333 

0.3333333 

5 

Probability of Capture 

9 

7 

1 

5 

5 

3 

9 

Reliability 

5 

5 

0.2 

1 

3 

3 

9 

Dynamic effects of capture and release 

3 

3 

0.2 

0.3333333 

1 

0.3333333 

9 

Minimal tether/payload considerations 

3 

3 

0.3333333 

0.3333333 

3 

1 

9 

Const. Cost 

0.1428571 

0.2 

0.1111111 

0.1111111 

0.1111111 

0.1111111 

1 


Sum 


Sum 


24.14 19.53 2.10 7.18 


12.78 8.11 49.00 


Minimal 

Dynamic tether/paylo 
effects of ad 


Current Probability capture and considerati 

Mass TRL of Capture Reliability release ons Const. Cost 


0.0414201 

0.0170648 

0.0529501 

0.0278638 

0.026087 

0.0410959 

0.1428571 

0.1242604 

0.0511945 

0.0680787 

0.0278638 

0.026087 

0.0410959 

0.1020408 

0.3727811 

0.3583618 

0.4765507 

0.6965944 

0.3913043 

0.369863 

0.1836735 

0.2071006 

0.2559727 

0.0953101 

0.1393189 

0.2347826 

0.369863 

0.1836735 

0.1242604 

0.1535836 

0.0953101 

0.0464396 

0.0782609 

0.0410959 

0.1836735 

0.1242604 

0.1535836 

0.1588502 

0.0464396 

0.2347826 

0.1232877 

0.1836735 

0.0059172 

0.0102389 

0.0529501 

0.0154799 

0.0086957 

0.0136986 

0.0204082 


1.00 1.00 1.00 1.00 1.00 1.00 1.00 


Ranking 

0.05 

0.06 

0.41 

0.21 

0.10 

0.15 

0.02 


Mass 

Current TRL 
Probability of Capture 
Reliability 

Dynamic effects of capture and release 
Minimal tether/payload considerations 
Const. Cost 


Code name bowfish bottleneck cablegrab boeing mini-v solid boom hayrake V-g ripper 




Mechanism 

.T08.P02a 


.T08.P02b2 


.T08.P02b1 


.T08.P02C 


T03d1.P04 


T03d2.P04.P 

T03d3.P04 


Criteria 

Weighting 

Score 

Raw Score 

Wtd Score 

Raw Score 

Wtd Score 

Raw Score 

Wtd Score 

Raw Score 

Wtd Score 

Raw Score 

Wtd Score 

Raw Score 

Wtd Score 

Raw Score 

Wtd Score 

Mass 

0.05 


8 

0.3992444 

8 

0.3992444 

8 

0.3992444 

7 

0.3493388 

3 

0.1497166 

2 

0.0998111 

4 

0.1996222 


















Current TRL 

0.06 


3 

0.1888376 

4 

0.2517834 

5 

0.3147293 

4 

0.2517834 

4 

0.2517834 

3 

0.1888376 

5 

0.3147293 


















Probability of Capture 

0.41 


3 

0.1888376 

5 

0.3147293 

6 

0.3776751 

5 

0.3147293 

4 

0.2517834 

4 

0.2517834 

5 

0.3147293 


















Reliability 

0.21 


4 

0.2517834 

3 

0.1888376 

4 

0.2517834 

2 

0.1258917 

3 

0.1888376 

1 

0.0629459 

4 

0.2517834 


















Dynamic effects of capture and release 

0.10 


5 

0.3147293 

5 

0.3147293 

7 

0.440621 

5 

0.3147293 

2 

0.1258917 

2 

0.1258917 

6 

0.3776751 


















Minimal tether/payload considerations 

0.15 


5 

0.3147293 

5 

0.3147293 

5 

0.3147293 

4 

0.2517834 

2 

0.1258917 

2 

0.1258917 

4 

0.2517834 


















Const. Cost 

0.02 


5 

0.3147293 

5 

0.3147293 

5 

0.3147293 

3 

0.1888376 

3 

0.1888376 

2 

0.1258917 

4 

0.2517834 


Total score (out of 10) 


1.9728908 2.0987825 2.4135118 1.7970935 1.2827421 0.9810531 1.9621062 Ta b[ e 

Table 4: Input Rankings for the AHP Analysis for Capture Mechanisms 
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Code name ugly stick Triangle Trap quad trap polysidedtrap scissor Iris bear trap and family zip tie original Iris 




Mechanism: 

T03d31.P04 

T074a.P03x 

T074b.P03x 

T074c.P03x 

T073 


T075 


T076 


T072 


Criteria 

Weighting 

Score 

Raw Score 

Wtd Score 

Raw Score 

Wtd Score 

Raw Score 

Wtd Score 

Raw Score 

Wtd Score 

Raw Score 

Wtd Score 

Raw Score 

Wtd Score 

Raw Score 

Wtd Score 

Raw Score 

Wtd Score 

Mass 

0.05 


5 

0.249528 

3 

0.149717 

3 

0.149717 

2 

0.099811 

2 

0.099811 

2 

0.099811 

4 

0.199622 

2 

0.099811 




















Current TRL 

0.06 


5 

0.314729 

7 

0.440621 

7 

0.440621 

7 

0.440621 

6 

0.377675 

7 

0.440621 

4 

0.251783 

8 

0.503567 




















Probability of Cap 

0.41 


5 

0.314729 

7 

0.440621 

8 

0.503567 

8 

0.503567 

7 

0.440621 

6 

0.377675 

7 

0.440621 

8 

0.503567 




















Reliability 

0.21 


6 

0.377675 

4 

0.251783 

4 

0.251783 

3 

0.188838 

4 

0.251783 

4 

0.251783 

2 

0.125892 

4 

0.251783 




















Dynamic effects < 

0.10 


6 

0.377675 

6 

0.377675 

6 

0.377675 

6 

0.377675 

6 

0.377675 

5 

0.314729 

3 

0.188838 

5 

0.314729 




















Minimal tether/pa 

0.15 


7 

0.440621 

5 

0.314729 

5 

0.314729 

5 

0.314729 

5 

0.314729 

4 

0.251783 

5 

0.314729 

4 

0.251783 




















Const. Cost 

0.02 


5 

0.314729 

3 

0.188838 

4 

0.251783 

3.5 

0.220311 

2 

0.125892 

4 

0.251783 

4 

0.251783 

2 

0.125892 


Total score (out of 10) 2.389687 2.163984 2.289876 2.145551 1.988187 1.988187 1.773269 2.051133 


Code name Balloon net turkey trap w/ harpoon Hamburger bun lacrosse Inflated net inflated torus Ball in Socket 




Mechanism: 

T03e.P05a 

T03.P05b 


T08.P05b 


T05c.Pxx 


T03e.P06a 

T06.P03 


T08.P06b 


Criteria 

Weighting 

Score 

Raw Score 

Wtd Score 

Raw Score 

Wtd Score 

Raw Score 

Wtd Score 

Raw Score 

Wtd Score 

Raw Score 

Wtd Score 

Raw Score 

Wtd Score 

Raw Score 

Wtd Score 

Mass 

0.05 


4 

0.199622 

6 

0.299433 

7 

0.349339 

3 

0.149717 

4 

0.199622 

8 

0.399244 

2 

0.099811 


















Current TRL 

0.06 


7 

0.440621 

6 

0.377675 

6 

0.377675 

4 

0.251783 

7 

0.440621 

2 

0.125892 

2 

0.125892 


















Probability of Cap 

0.41 


4 

0.251783 

3 

0.188838 

4 

0.251783 

4 

0.251783 

4 

0.251783 

2 

0.125892 

2 

0.125892 


















Reliability 

0.21 


7 

0.440621 

4 

0.251783 

4 

0.251783 

3 

0.188838 

7 

0.440621 

2 

0.125892 

3 

0.188838 


















Dynamic effects < 

0.10 


6 

0.377675 

6 

0.377675 

6 

0.377675 

2 

0.125892 

6 

0.377675 

3 

0.188838 

3 

0.188838 


















Minimal tether/pa 

0.15 


4 

0.251783 

4 

0.251783 

4 

0.251783 

4 

0.251783 

4 

0.251783 

4 

0.251783 

3 

0.188838 


















Const. Cost 

0.02 


4 

0.251783 

7 

0.440621 

7 

0.440621 

4 

0.251783 

4 

0.251783 

4 

0.251783 

3 

0.188838 


Total score (out of 10) 2.21389 2.187809 2.30066 1.47158 2.21389 1.469324 1.106945 

Table 5: Weightings and Scores for the AHP Analysis for Capture Mechanisms 
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4.0 BASIC DESIGN SELECTION AND DEVELOPMENT 

(For Ideal Capture Device) 

Capture Mechanism Design 

The brainstorming sessions gave birth to a particularly large number of capture mechanisms deemed 
worthy to perform a more in depth study of. A first order trade study was performed on these 
mechanisms and the four most promising designs were chosen to endure a brief two week engineering 
evaluation. These four capture mechanisms were to be presented to the evaluation team of Lockheed 
Martin/NASA. Of these four, the top mechanism was to be thoroughly designed and a prototype 
assembled, tested, and evaluated. This process was completed with the prime candidates for further 
development being: 

1 . Modified Iris (quad trap) 

2. Turkey Trap (Seth noose) 

3. Tethered Flyer 

4. V-Gripper 

This section will provide a first order engineering analysis of each of these four designs (contained in 
Section A 4.1 through 4.4). These sections immediately follow. The results and conclusions of a 
study of these first order analyses are contained in Section A 4.5. Based on the conclusions from 
Section A 4.5, a more detailed analysis of the Quad trap mechanism is conducted in Section A 4.6. 
Finally, an analysis of an additional mechanism for MXER tether systems, a 3 dof mechanism for 
pointing solar arrays is provided in Section A 4.7. 

4.1 Modified Iris(Quad Trap): First Round Design Analysis 

4.1.1 Overview 

This section will summarize the results of a first-round of design and analysis of the modified iris 
capture concept. The primary elements of the design are described, as well as the analysis procedures. 
The key results are summarized and conclusions are made regarding suitability and future development 
of this capture concept. 

4.1.2 Brief Description 

The quad trap is designed to accept an extension from the payload and direct it to the corner plate. For 
this analysis, that extension is a single cable with a gyroscope at one end for stabilization of the 
trajectory. This lanyard can be reeled in and reloaded to use as many times as needed, using only 
mechanical energy for propulsion. The quad trap has smart sensors on it to detect the presence of the 
lanyard and, upon entry of the lanyard into the capture area of the quad trap, the trap actuates and the 
lanyard is forced into a plate located at one of its four corners. This plate is designed to hold the full 
force of the payload as it is accelerated by the tether. 

4.1.3 Summary of Analysis Performed 

Several key analyses were perfonned to characterize system behavior and analyze results. These were: 

1. Time closure for the iris upon entry of the lanyard into capture area were calculated to be around 
three and a half seconds for gravity assisted closing. 

2. Stress calculations on spreader bars determined. 

3. Sizing of the spreader bars(0.67 slugs) and comer plates(1.86 slugs) was perfonned. 

4. Vibration calculations to detennine natural frequency of spreader bars were perfonned to 
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ensure the tether was not driving the capture system.(f n =662 kHz) 

5. Design for reset perfonned and suggested ratchet designed. 

6. Mass for lanyard gyroscope determined. 

4.1.4 Summary of Key Conclusions 

1 . The closure of the system can be perfonned entirely by inertia keeping the system simple. 

2. The parameters of the system suggest high probability of capture and part reliability. 

3. A tethered lanyard was chosen to connect the payload to the iris mechanism due to its 
simplicity and feasibility. This lanyard could be light and easily packaged on the payload. 

4. The capture plate has been designed to receive the full load of the payload co-aligning the 
stresses with the cables making the configuration attractive. 

5. The design of this capture device should be pursued based on simplicity of design and 
reliability of mechanism both in capture and in performance. 


As noted, this mechanism is capable of closing within the allotted closure time frame entirely from 
inertial effects, making the closing process entirely passive. This can be demonstrated by assuming a 
simple pendulum model. Neglecting frictional effects, the period for a pendulum is given by the 
following relation, 



d(f> 

l-sin 2 (— )sin 2 (^) 



where K denotes the elliptic integral and the factor 2KJ n < 1 .002 for 0 < 10° . Tabulated values for 
elliptic integrals are available in standard references, for example. With a value of L of 100 meters to 
the end of the tether and g replaced by the local acceleration of 19.8 m/s , the period x is found to be 
14.12 seconds. One-quarter of this period is the iris closure time, about 3.53 seconds. Preliminary 
simulations suggest that the time interval during which the payload is within the capture window is 
about six seconds so that the inertially-driven closure rate would suggest full closure perhaps one-half 
second after CAP. 


For the release stage, the payload will jettison its lanyard severing the connection between payload and 
tether. The remaining lanyard will then be shed from the Marshall plate using the lanyard shed 
mechanism which will push the payload lanyard away from the groove in the Marshall plate. Once 
shed, the payload lanyard will fall away from the tether due to centripetal acceleration. In this way, 
the trajectory of the spent payload lanyard can be controlled. Figure 4. 1.4-1 below is a conceptual 
drawing of the quad trap. 
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Tether attachment points 


Open/closed 

configurations 



Motor-driven 
opening jacks 


Corner plate Spreader bars 

Mass=1.86 slugs Mass=0.67 slugs 

Material: titanium-alloy Material: titanium-alloy 

Inertial acceleration (g): 2x32.2 

Figure 4.1.4-1: Schematic of the Iris (Quad Trap) Mechanism Concept 



T.»,,=ing/cose 



EM = — [I]m 
dt 


EF = — pv = ma 
dt 


Sum of forces and moments with simplified calculations for scaling magnitude 
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(a) 

Open configuration of iris 


(b) 


One proposed cross section of 
spreader bar 

Figure 4.1.4-2: Open Configuration of the Iris and Cross Section of the Spreader Bar 


As Tbar is the reaction needed to maintain the “open” configuration, it is calculated by a summation of 
forces on the block and is found to be mg tan(O). Assuming the corner blocks each have an earth 
weight of 60 lbs and the “spreader bars” have an earth weight of 80 lbs each, the maximum load Tb ar 
would be 4.8 lbs, which is the axial load carried by the “spreader bars.” Based on Euler buckling 
theory, a critical buckling load for pinned ends is given by, 

n - EIk1 
“ L 2 ’ 

where F cr is the stress that will cause the member to buckle, E is modulus of elasticity, I is the moment 
of inertia, and L is the length of the beam. Assuming that this member has the cross section shown and 
that the length is fixed, one can see that F cr is driven by E and 7. The material would likely be an I- 
beam cross section of some titanium alloy. 


I={bd 3 -h 3 (b-t)}/12 

Using material properties for titanium (E approximately 15 times 10 6 psi), an area moment of inertia of 
7=1 7in 4 and a length L of 394 in, the critical buckling load is approximately 16 kips, giving a safety 
factor of more than 3000. Additionally, the spread bars should not be expected to exhibit low 
frequency vibrations. For pinned ends, frequencies of vibration are given by, 



so that the fundamental frequency is 662 kHz. 


The capture plate Figure 4. 1 .4-3 rests on a corner of the quad trap and takes the full load once the 
satellite and the tether connect. The maximum loads to be carried should be approximately 10,0001bs. 
When this static load was applied to the plate in an FEM 
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model, the maximum stresses in the plate were found to be approximately 10 ksi, well under the yield 
stress for some forms of titanium-alloy (130-140 ksi). The large stress margin provides significant 
reliability even though the actual loading will be dynamic. 



Figure 4.1.4-3: FEA of Results for Maximum Loads on the Capture Plate 

After the quad trap closes due to inertial forces, it resets itself to the open configuration using a type of 
ratchet system. Within the ratchet are two arms, one for resetting the spreader bar to the open 
configuration and one for holding it in place. Small bearings are mounted inside the casing for this 
ratchet system to aid in closing. Upon arrival of the payload lanyard within the designated capture 
area, a small servo behind the holding ratchet pushes it forward releasing the bar to slide through. The 
reset ratchet is run by a small motor which only needs to provide about .4 ft-lbs of torque. The 
Payload Lanyard will consist of a small tethered gyroscope cast out from the payload satellite using a 
preloaded spring. 

Prior to the time for deployment a small gyroscope attached to the end of a -420m satellite mounted 
tether will be spun up. At time for deployment this gyroscope will be launched from the satellite 
toward the tether via a pre-loaded spring. As the gyroscope flies it will pull the lanyard from its spool. 
This is demonstrated conceptually in Figure 4. 1.4-4 below. 
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Figure 4.1.4-4: Concept of the Gyro-Stabilized Payload Lanyard 

The gyroscopic effects of the spinning gyroscope will help to stabilize the trajectory of the lanyard as 
it deploys. Some number of stops will be woven into the top portion of the lanyard to allow it to 
engage the Marshall Plate. These stops will be wound onto the spool along with the rest of the lanyard 
and will deploy in the same fashion. When the lanyard reaches the limit of its deployment, the storage 
spool will provide braking torque to slow the gyroscope and prevent “bounce back” effects. 

At full deployment the gyroscopic affects of the end-mounted gyroscope will resist moments imposed 
on the lanyard by outside forces. For the payload lanyard assume: 

M out = ]n 3 

UJ R =1100 RPM a, =4 1469 12 rad/s 
M = 25g/m = 0.025kg/m 

otmax=0. 01 rad/s n 3 

(Lanyard approximates thin bar inertial properties) So the fonnulation for the inertia of the gyroscope 
proceeds as: 

I r =— = 8333.3 kg/m 2 
3 

I Gb 2 = 1 /2mR 2 = 1 /2m(. 1 0 1 6/2) 2 =0 .00 1 29m 

, -i-i 1016 i 

I bl =I b3 =\/4mR 2 +\/12mL 2 +md 2 =l/4m (— — ) 2 +1/12 m( 100) A 2= 10,000m 

cos(#) sin/#) 0 cos(^) 0 sin(^) 

R a n = -sin(#) cos(<9) 0 and R° = 0 10 

0 0 1 -sin(^) 0 cos(^) 
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then the lanyard inertia I z =R c a 8333a_3 = 8333 n3 

EM = la + 

and can readily be shown that 

_ -82.33 + 1.43M7sin(fl)cos(fl) 

1 0 + 1 .IE 1 7 sin(6*) 008(6*) 

These gyroscopic effects in combination with the tension provided by the braking reel will help to 
preserve the orientation of the lanyard for as long a period as possible. However the tension in the 
lanyard that aids with its orientation also causes the lanyard to move back toward the payload. There- 
fore as the lanyard slowly retracts back toward the payload the braking reel will be used to re-spool the 
excess slack. Once they lanyard is fully retracted back onto the spool the spool torque will be used to 
reload the launcher spring for another sequence. 

4.1.5 Other Elements for the Iris Mechanism under Consideration 
Shape Memory Alloy 

A strand of shape memory alloy is woven into the payload lanyard in such a way that when the lanyard 
is un-spooled the shape memory alloy is activated causing the lanyard to straighten and stiffen. This 
concept was pushed aside due to a general ignorance of the specifics of shape memory alloys. In the 
first case it is believed that the extreme temperature fluctuations may prevent the SMA from stiffening 
properly. It is also believed that a structure of sufficient stiffness would likely be quite heavy. 

Claw-Hammer Hinge 

The claw hammer hinge is a concept that would potentially replace the Marshall Plate as a means of 
making the terminal connection to the payload boom (Figure 4. 1.5-1). In the capture configuration the 
claw hammer hinge would be folded out ward in such a way that the stops in the payload boom would 
slide into the “claw hammer” shaped groove cut in the device. At the time for payload release the claw 
hammer hinge would then pitch backward allowing the payload boom to slide all the way through the 
claw hammer groove and release from the tether. 



Figure 4.1.5-1: Schematic of the Claw-Hammer Concept 
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4.2 Turkey Trap: First Round Design Analysis 

4.2.1 Overview 

This section will summarize the results of a first-round of design and analysis of the Turkey Trap 
capture concept. The primary elements of the design are described, as well as the analysis procedures. 
The key results are summarized and conclusions are made regarding suitability and future development 
of this capture concept. 

4.2.2 Brief Description 

The Turkey trap is a cable system, enclosing an area of 10m radius by 100m length (Figure 4. 2. 3-1), 
which attaches to the payload and extends out to meet a specially designed grapple on the end of the 
tether. The grapple for this concept is called the “Hamburger Bun” (Figure 4. 2. 3-2), an oblong 
spherical mechanism with retracting fingers attached to the end of the tether. Weights on the cables of 
the Turkey trap are ejected from the payload to the capture area and the mouth of the Turkey trap is 
expanded by inflatable pillars. After the Hamburger Bun passes thru the mouth of the trap, capture 
fingers within it are extended perpendicularly outward to meet the Turkey trap. Once this has occurred 
the inflatable’s are deflated and the noose constricts while being reeled toward the payload where it 
attaches to the Hamburger Bun. 

The concept behind Turkey trap is simple in that it contains a net and a type of hook which, when in 
contact, simply hook the net with a special type of prong. Upon closer inspection, it became apparent 
that solving the dynamics for the system was no trivial task. Although nets have been used for 
thousands of years, no governing equations exist, nor are likely to exist, for a body with such a large 
number of degrees of freedom. Eventually the design and tests moved from derivable to experimental. 
The evolution of the net was an area where some design could be applied, beginning with the standard 
cris-crossing of strands of cable as in Figure 4.2.3- 1 and the grappling device in Figure 4.2. 3-2 below. 


(a) (b) 

Original turkey trap Modified "laundry bag" net with drawstring motors 

Figure 4.2.3-1: Schematic of the Turky Trap Concept 
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(a) (b) 

Hamburger capture device Model of internal workings of hamburger device 

Figure 4.2.3-2: Schematic of the Hamburger Bun Grapple Concept 

Upon inspection of this original concept, it was noted that the probability of capture was highly chaotic 
in that not all hooks might catch a portion of the net. If such a scenario occurred, the worst case 
scenario would provide one prong hooking on one strand. This would dictate that every strand would 
need to be robust enough to carry the entire weight of the accelerated load. This design was 
unacceptable as the total weight of the net would be too much, and the transportation of it too difficult. 
To correct for this inefficiency, a drawstring system of sorts could be implemented to heighten the 
probability of capture. Some number of nodes were arranged around the periphery of the mouth of the 
net (Figure 4. 2. 3-3). Upon entry of the capture device into the net, the drawstrings would close rapidly 
ensuring that all prongs on the capture device would interact with the netting. Originally, the net was 
assumed to be the full volume of the error window, making it 10 meters by 100 meters. To decrease 
the weight of the net, the size of the net was diminished to simply the mouth and a small extension as 
in Figure 4. 2. 3-3 below. 



(a) (b) 

Truncated net with drawstring motors Noose-type Turkey trap 

Figure 4.2.3-3: Modified Versions of the Turkey Trap 


This net would also have the drawstrings on the mouth, but would also have some type of reel at the 
base. When the capture concept entered into the mouth, the drawstrings would close upon the device. 
Any overshoot on part of the tether tip would be countered by the reel at the base of the satellite. The 
next generation of concepts was radically different in that the net was broken down into its basics; it 
needed to be lightweight and it needed to catch at a known number of points. To achieve this, the 


31 



Noose-type Turkey trap was suggested. Originally it would consist of simple sliders on a number of 
stings. The net system would begin deployment by spinning the net. The sliders, having some mass 
would be ejected in the direction of the capture location still spinning, causing them to move away 
from each other widening the mouth of the Noose-type Turkey trap. The cables would go from the 
base of the payload out to one slider, loop around to another slider, and back to the payload. When the 
capture device went thru the mouth of the noose, the cables would be retracted suddenly. Momentum 
would keep the sliders from speeding up as fast as the cables, and consequently, the sliders would 
follow the path around the top of the cable arc traveling toward each other. A model of this was built 
and tested, and found to be somewhat difficult to predict. 

4.2.3 Summary of Analysis Performed 

After briefly engaging in modeling the net, it was concluded that the system could not be modeled 
practically due to the large number of degrees of freedom associated with the net. Based on the nature 
of the system, the suggestion was made to move to an experimental simulation approach. An 
experimental model was built with the intent to discover any correlations between geometric 
parameters and system behavior. Further details on the experimental analysis are provided in Section 
4.2.5 below. 

4.2.4 Summary of Key Conclusions 

1) The Turkey trap was found to be an attractive net as it minimizes weight while increasing the 
probability of capture. 

2) Several inflatables would form a perimeter of the Noose-type Turkey trap for entry of the 
Hamburger Bun. A minimum compressive force would be needed to hold an open configuration. 

3) Limited capture attempts exist in the current configuration due to inflation requirements. 

4) This concept has great potential but also significant modeling costs due to its reliance on nets and 
inflatables. 

During the assigned two week design study phase the generalized turkey trap concept evolved thru 
three different configurations. The first such configuration has been referred to as the “laundry bag” 
the second configuration “the Noose” and the third configuration “the Sethnoose”. All of these 
configurations are intended to mate to some variation of the “hamburger Bun” mace concept. 

4.2.5 Experimental Analysis 

Due to the non-linear nature of the equations of motion as applied to the Noose-type Turkey trap, any 
model must be quite complex before it can be considered to be representative of the system. In order 
to achieve some representative results in the required time span an experimental approach was 
adopted. This experimental approach would investigate the stability and controllability of the proposed 
system. 

The experiments were conducted as follows. First, a model was constructed. This model consisted of 
7 sections of aluminum rod which were drilled with holes so that they were somewhat representative 
of the pulley masses proposed in the Noose-type Turkey trap system. These pulley masses where then 
strung together using lightweight fishing line in a manner similar to that proposed by this system. The 
experimental model was then set spinning using a variety of different techniques to test its stability and 
controllability. First the model was spun using a simple cordless drill attached to the point where the 
lines all converge. This particular setup produced mixed results. Initially the pulley masses would 
align into a circular configuration although they did not extend radially as far as was hoped. After a 
period of about 15 seconds the circular configuration would then deteriorate into a second stable, but 
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non geometric configuration where it would remain until stopped. These results were remarkably 
repeatable from test to test and showed a repeatable time-to-decay. 

Several hypotheses were proposed to explain the mixed results of this first experiment. Most 
importantly the experiment occurred in a gravitational system. In a gravitational system the pulley 
masses would need a much higher velocity to extend into a large radius circle because the centrifugal 
force has to overcome gravity. A second hypothesis centered on the irregularity in friction between 
the fishing lines and the pulley masses. Given that each pulley mass is only a drilled hole the friction 
thru that hole is assumed to be quite irregular from part to part. Once the system gains enough speed 
for the pulley masses to begin sliding apart from one another any unbalance in friction between one 
pulley mass and the next could cause instability. Pulley masses were initially not spread apart from 
one another and it was proposed that the initial condition of no radius resulted in an instability. There 
were also concerns that the means of driving the system caused a substantial amount of twist in the 
lines which may have caused the instability due to a difference in line length. 

The next experiment of note involved rotating the system from a turntable. This experiment was 
intended to address the concerns of line twisting and zero initial circle radius. Pulley masses were 
arranged on a turntable in a circular configuration some distance from one another. The turntable was 
then rotated to impart a spinning speed to the system and the pulley masses were lifted off of the 
turntable. This experiment was somewhat more successful. If the speed of the turntable was well 
controlled the system would fall into a stable circular configuration and remain in that circular 
configuration for some time. However, abrupt changes in rotation speed would disturb the system and 
cause it to decay into a random configuration. The same hypotheses on failure listed above also apply 
here. 

The final experiment used drinking straws to test the hypothesis that some residual stiffness in the lines 
would dampen instabilities in the system and provide an initial radius. Drinking straws were applied 
to the lines of the Noose-type Turkey trap so that they divided the pulley masses by a distance of ~3”. 
The model was then spun up using the same cordless drill arrangement as described in the first test. 
Results were much better in this final test. Although the pulley masses did not spread apart much 
farther than the straw sections had forced them to the configuration was much more stable and could 
be sustained for quite some time. 

It was then suggested that, instead of a rotation to keep the weights in place, an inflatable was used as 
shown in Figure 4.2.4- 1. This would ensure that the sliders were in an appropriate configuration when 
the capture device entered. Upon entry of the capture device, the inflatables would deflate, the cables 
tugged back toward the satellite, and the sliders move to the center to meet the capture device. 
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Figure 4.2.4-1: Noose-Type Turkey Trap with Inflatables 
4.3 Goalpost Flier: First-Round Design Analysis 

4.3.1 Overview 

This document will summarize the results of a first-round of design and analysis of the goalpost flier 
capture concept. The primary elements of the design are described, as well as the analysis procedures. 
The key results are summarized and conclusions are made regarding suitability and future development 
of this capture concept. 

4.3.2 Brief Description 

The goalpost flier capture concept, shown in Figure 1, consists of a “goalpost” on the tether (a set of 
light booms that guide the cable toward a cable grab) and a small tethered passive flier on the payload. 
The tether flier is launched on a simple trajectory and places the cable over the goalpost and into the 
cable grab. The flier will maintain tension in the cable until the cable grab pulls taut (occurs when the 
relative velocity b/n the tether and payload goes to zero). The tethered flier is then cut free and travels 
harmlessly into the confines of deep space. See Figure 4.3.2- 1 below. 



Figure 4.3 .2-1: Goal-Post Filter Schematic 

4.3.3 Summary of Design Points 

The goalpost on the tether consists of composite booms that terminate in the base of the cable grab as 
shown in figure 2. The base of the tether is encased in a shield to protect it from impact with the 
tethered flier. The tethered flier is passively launched and can fly a fixed trajectory or a variable 
trajectory (if given a small thruster). The flier-cable consists of 20 mm HMPE wire rope. The 
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grasping mechanism is a cable grab, a cam device that allows the cable to travel in one direction only. 

An additional tether shearing mechanism is added to release the tethered flier after grasp. 

4.3.4 Summary of Analysis Performed 

Several key analyses were perfonned to characterize system behavior and analyze results. These were: 

1 ) Dynamic simulation of the interaction between the flying cable and the light boom “goal 
posts”. This simulation assumed a small tension in the flying cable (approx. 1 N) and was 
performed in ADAMS using thin beams and line elements. The simulation demonstrated the 
desired performance, the flying wire would be guided toward and into the cable grab with 
minimal tension. Samples from the dynamic simulation are shown at the end of this report. 

2) Size of the flying wire are determined. These are given as HMPE 20 mm dia. wire rope, total 
strength of 1 17.5 kN, weight of . 188 kg/m. 

3) Mass estimate of the system was performed. The flier is estimated at 30 kg, the flying tether at 
30 kg, the goalpost booms and cable grabs at 190 kg. Total mass estimate is 250 kg. 

4.3.5 Summary of Key Conclusions 

In light of the resulting analyses, several conclusions are drawn. 

1 ) From a dynamic standpoint, a small amount of tether tension will guide the tether along the 
goal post and into the cable grab. Other dynamic effects of impact do not seem to be limiting. 

2) An appropriate trajectory needs to be detennined. Currently, simple rules in defining this 
trajectory can be constructed. A curvilinear trajectory (an accelerating trajectory) is most 
likely necessary. Therefore, a completely passive flier may not be possible. 

3) The tether itself will always exist in the field of one or more potential trajectories. It is 
proposed to shield the end of the tether and allow the possibility of collision b/n the tether and 
tethered flier. More modeling will be performed to demonstrate that this is not a problem. 

4) The mass of the system is low, the probability of the system working seem reasonable. The 
complexity of the system is low to moderate. 

5) Plans for release of the payload and release of the tethered flier are still needed. 

6) Due to the low mass and low to moderate complexity, this design should be pursued further. 

Figures 4.3.5- 1 through 4. 3. 5-4 show a simulation of the tethered flier (performed in ADAMS) making 

contact with the goal post and moving toward the cable grab. 
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Figure 4.3.5-1: 


Figure 4.3. 5-2: 
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Simulation Results of Contact between Tethered Flier and Goal Posts, t=0.0 



Simulation Results of Contact between Tethered Flier and Goal Posts, t=0.93 
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Figure 4.3.5-3: Simulation Results of Contact between Tethered Flier and Goal Posts, t=1.4 



Figure 4.3.5-4: Simulation Results of Contact between Tethered Flier and Goal Posts, t=3.7 
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4.4 Modified V-Gripper: First-round Design Analysis 

4.4.1 Overview 

This document will summarize the results of a first-round of design and analysis of the modified V- 
Gripper capture concept. The primary elements of the design are described, as well as the analysis 
procedures. The key results are summarized and conclusions are made regarding suitability and future 
development of this capture concept. 

4.4.2 Brief Description 

The modified V-Gripper, shown in Figure 4.4.2- 1 below, consists of a planar, guiding mechanism on 
the tether and an extended trap wire on the payload. The guiding and gripping mechanism on the 
tether spans an area approximately 20m x 100m, and is actuated by a cable in the center of the 
mechanism. As the payload and tether approach rendezvous, the V-gripper retracts to guide the trap 
wire toward the center of the V-gripper, and to complete the grip when the payload and tether have 
zero-relative velocity. The V-gripper has zero-degrees of freedom. Rather it relies on elastic 
deformation to form the guiding surface and to retract as the trap wire approaches the tether. 



Figure 4.4.2-1: V-Gripper Schematic 

4.4.3 Summary of Design Points 

The modified V-gripper consists of a composite rod with annular cross-section shown in figure 2. The 
trap wire supported on the payload consists of 20 mm HMPE wire rope and two light-weight booms. 
The cable grasping mechanism is a single dof, two jaw grasping device, made of composite material. 
The V-gripper is driven by a 20 mm HMPE wire rope, approximately 150m in length with actuated 
cable drum. 

4.4.4 Summary of Analysis Performed 

Several key analyses were performed to characterize system behavior and analyze results. These were: 
1. Dynamic simulation of trap wire impacting and sliding along the V-gripper. Simulation was 
performed in ADAMS using thin-beam sections to form the model. The simulation 
demonstrated the desired performance of the trap wire guided toward the center of the gripper 
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with very minimal contact force (sub 1 N force). Screen grabs of the dynamic simulation are 
shown at the end of this report. 

2. Strain analysis was performed on the V-Gripper during a full deflection. Max strains in the 
device were very small (sub 1%). Some of the results of this analysis are shown at the end of 
this report. 

3. Deflection analysis was performed on V-Gripper during normal operation. Nonnal operation 
places the device in an approximately 2g loading condition. Due to the light-weight nature of 
the beam, the 2g load causes large deflection of the device. Off-loading cables will be 
necessary to allow V-gripper to function. 

4. Size of the trap-wire and actuating wire are detennined. These are given as HMPE 20 lmn dia. 
Wire rope, total strength of 1 17.5 kN, weight of . 188 kg/m. 

5. Mass estimate of the system was performed. The primary contributors to the mass were: 1) the 
V-gripper boom (approx. 250 m in length) 2) HMPE cable. The total mass of this device in an 
ideal case appears to exceed 500 kg and perhaps 1000 kg. 

4.4.5 Summary of Key Conclusions 

In light of the resulting analyses, four conclusions are drawn. 

1. From a dynamic and operational standpoint, the modified V-Gripper device is demonstrated to 
work as expected. Impact b/n trap wire and V-gripper Boom does not create instabilities, 
dynamic response shows a clear guiding of the trap wire toward the center. 

2. The strain in the boom is minimal to deflect to the desired shape. 

3. Due to the accelerating state of the end of the tether, approx. 2g, the V-gripper boom is under 
significant quasi-static load. This can be accounted for with off-loading wires (analysis is 
performed to demonstrate this can work). However, the trade-off is added mass and 
complexity. 

4. The mass of the system is exceedingly high. The long boom, off-loading wires, and long trap 
wire contribute to this high mass budget. 

5. Due to the high mass estimate and additional complexities of off-loading, this concept will not 
be further pursued for the current time. 

4.4.6 Finite Element Analysis of Concept Name: V-GRIPPER 

A simple model of the V-Gripper is considered in a coimnercial finite element analysis software 
(ANSYS). The booms are 4 in diameter and made of graphite. Material specifications are: E= 1 0 Gpa 
and p=1800 kg/m with dimensions shown in Figure 4.4. 6-1. Following this, a series of figures depict 
typical responses of the system for various loads or displacements (Figures 4.4. 6-2 through 4.4.6-13 
below). 
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Figure 4.4.6-1: Sketch of Vi of the Ugly Stick V-Gripper 



Figure 4.4.6-2: Displacement of Half of V-Gripper 

(With 1 N load applied as shown) 
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Figure 4.4.6-3: Displacement Profile 

(With an input displacement of 100 m at the center) 



Figure 4.4.6-4: Strain Distribution for Half of V-Gripper 

(With 100 m input displacement) 
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Figure 4.4.6-5: Stress Distribution for Half of V-Gripper 

(With 100 m input displacement) 
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Figure 4.4.6-6: Displacement Plot 

(Less than 1 00 m input for full case) 
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Figure 4.4.6-7: Strain Distribution 

(For 100 m displacement at center) 



Figure 4.4.6-8: Stress Distribution 

(For 100 m displacement at center) 



Figure 4.4.6-9: Displacement with 1 N load 

(Applied off center as shown) 
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Figure 4.4.6-10: Stress Distribution 

(For Figure 4. 4. 6-9) 
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Figure 4.4.6-11: Deflection with 1 N Load 

(Applied as shown) 
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Figure 4.4.6-12: Stress Distribution 

(For Figure 4.4.6- 1 1) 
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Figure 4.4.6-13: Deflection for Case of 1 N Load in Center 

4.5 Trade Study on Second Round Selected Designs 

After a brief study of tradeoffs on the top four choices, the Modified Iris was selected as being the 
simplest, most capable idea under the conditions and stipulations placed on the design process. One of 
the most attractive features being the passive gravity-assist closing of the device. It was noted that 
under the force of its own weight, the mechanism could close within an appropriate time frame. This 
was a colossal advantage over the active mechanisms or mechanisms requiring more complex methods 
or aid in closing. Upon a firm agreement and support of team members and approval from the project 
manager the modified Iris (Quad Trap) was designated the candidate which would be subjected to 
physical testing. 

4.6 Develop Advanced Models of Selected Mechanism (Quad Trap) 

Upon choosing an appropriate candidate for further investigation, a second order engineering analysis 
was performed upon the more crucial elements of the system. The exact geometry was defined and an 
ANSYS model was created to observe stresses and deflection due to loading. Over several weeks, this 
model was shaped and manipulated to meet weight/strength criteria with several tests performed 
permutations on shape and material. The final material choice was an intuitive guess, a titanium alloy, 
and the final shape can be seen in Figure 4.6-1. 



Figure 4.6-1: View of Catch Plate 
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The evolution of the plate was drawn from the criteria set up declaring function and weight 
requirements. The original concept was a method for moving weight from off axis of the tether onto a 
collinear position. The plate concept was the initial choice, followed closely by a set of hooks, see 
Figure 4.6-2 below. 



Figure 4.6-2: Capture Hooks to Replace Plate 

The hooks were cast aside from the uncooperative nature of the geometry. The initial plate is shown in 
Figure 4.6-3 which handled stresses appropriately but held an expensive weight penalty. 



Figure 4.6-3: Initial Plate Response to Typical Loads 

The initial plate handled the an assumed load often thousand pounds with resulting stresses at lOksi, 
well under the limit of yield stress of titanium, 25ksi. The redesign of the plate brought its weight to 
lOOlbs with a max stress of 12ksi which, while a significant drop in weight from the original 2001bs, 
was still a non-trivial weight penalty and high stresses as well. The intermediate designs are shown in 
Figure 4.6-4 below. 



(a) 


(b) 
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(C) 

Figure 4.6-4: Intermediate Degins of Capture Plate 

The final design of the plate brought its typical stresses down to under lOksi with maximum stresses at 
13ksi, but with the weight peaking at 601bs. This was well within the design criteria and a suitable 
recommended design. The next design of the Modified Iris (quad trap) was to attach the plates to 
sliding railings. The railings were designed as I-beams for their strong moment resistance while 
maintaining very low weight. A picture of the final design can be seen in Figure 4.6-5. 



Figure 4.6-5: Final Second-Order Design of the Modified Iris 

One of the other mechanisms which was addressed was the payload side of the iris. It was suggested 
that the payload side of the capture mechanism be simplified from the extendable boom to a single 
cable which would have nodes on it to be grasped by the capture plate on the modified Iris. A 
casting/reeling mechanism was designed to launch the lanyard to be captured by the quad trap. One of 
the desired features of the launcher was that it must have no limit to the number of launches, deeming 
all forms of propellant out of the design race. The launcher could also be packaged in a format to 
which satellite manufacturers could simply attach it to their satellite pre-launch with no hassles. See 
Figure 4.6-6 below. 
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Figure 4.6-6: Schematic of Lanyard Launch Mechanism 

Several models were drawn up in ADAMS dynamic modeling software during the process as well. A 
simplified version of the quad trap was created to aid in the design of the beams. The beam was 
initially modeled after an Able boom, a technology which had been discussed a great deal and 
generally accepted. From this model the maximum deflection of the central guidance beam was 
determined. The model consisted of two blocks (both steel 2’ x 2’ x 2”) and a circular 10m long beam 
(8’ diameter) modeled after one of the Able Booms. See Figure 4.6-7 below. 



Figure 4.6-7: Simplified Version of Deloyable Boom 


The two blocks were simulated free to move along the pole but constrained from moving off the pole. 
A force of 5600N , a force assumed to occur under 2G loading, was placed on the center of mass of 
each of the blocks forcing them closed as well as another 5600N pushing the entire system 5m 
horizontally due to the symmetry of the system. This horizontal motion caused the beam to deflect and 
the maximum deflection with respect to time was calculated parametrically with variable beam 
stiffness and mass, and three different cases simulated for each. The forces required to close the 
mechanism in one second was also calculated for each of the masses. Figure 4.6-8 summarizes the 
ADAMS model of the Quadtrap system. 
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Figure 4.6-8: Schematic of ADAMS Model of Quad Trap 

A summary of the results are shown in Tables 6-9. These tables demonstrate the system response in 
terms of displacement, velocity, system energy and acceleration respectively. 


Table 6: Displacement of Quad Trap during Capture 



Analysis: Last_Run 


Time (sec) 


2004 - 02 - 2414 : 43:34 
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Velocity (inch/sec) 


Table 7: Velocity of Quad Trap During Capture 


model 1 



Analysis: Last_Run Time (sec) 2004-02-2414:43:34 


Table 8: Sum of System Energy in Quad trap 


model 1 
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Table 9: Acceleration of Quad Trap During Capture 


model 1 



After these intermediate models were created, the scaled quad trap prototype was modeled. The 
developed prototype is discussed in the section A 5. 

4.7 Evaluation of 3 DOF Spatial Joint Mechanism 

(For MXER and other in-space propulsion concepts) 

4.7.1 Summary of the 3D Solar- Array Tracking Joint Design for Solar- Array Tracking 

A three degree-of-freedom, spatial, parallel-architecture manipulator (referred to as a spatial joint) is 
investigated here as an improved gimbal mechanism for in-space applications. This spatial joint has 
many advantages, including large, greater-than hemispherical, singularity-free workspace, ground- 
relative actuation, and a high strength-to-weight ratio. For this effort, two prototype spatial joints were 
designed, fabricated and delivered. These prototypes demonstrated real-time closed-loop control by 
tracking an operator input via ajoystick and are shown in Figure 4.7. 1-1 (a) and (b) below. 



(a) (b) 

Figure 4.7.1-1: Prototype of 3 DOF Mechanism in Two Positions 


52 




As an example of implementation, the design of this spatial joint for implementation as the primary 
gimbal support for a sun-tracking solar array on the CEV. The requirements for this joint are to 820 ft- 
lb torque load while providing tracking capability to the solar arrays (wrist pitched to 90 degrees while 
rotating fully). Figures 4.7. 1-2 and 4.7. 1-3 provide an overview of the design as a CAD model. 



Figure 4.7.1-2: CAD Rendering of Solar Arrary Tracking Joint 



Figure 4.7.1-3: Solar Array Tracking Joint 

(At 90 Degree pitch at location in full rotation) 

In Figures 4.7. 1-2 and 4.7. 1-3 above, the Solar- Array is removed for clarity. The location in Figure 
4.7. 1-3 above shows the joint in the highest-load configuration. The actuation system is currently 
demonstrated as a combination Maxon EC-90 Flat motor and HDC-1M harmonic drive gear head. 
Numerous configurations for actuation are available; this one is shown for representation only. 
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Figure 4.7. 1-4 demonstrates the overall kinematic parameters for the current joint design to achieve a 
greater-than hemispherical workspace. With a firmer set of requirements and more detailed design, 
this workspace can be scaled to a reasonable degree. 
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Figure 4.7.1-4: Wire Frame Solar Array Tracking Joint 

(With overall dimensions) 


Key metrics for the initial design of the 3D solar-array tracking mechanism are listed in the following 
table: 


Table 10: Summary of 3D Solar-Array Tracking Mechanism 


* 


wit 


Component 

Description 

Mass/ Workspace 

Motor Drive: 

Maxon EC90 Flat Motor 

1.5 lb 

Drive reduction 

Harmonic Drive HDC 1M 

2.61b 

Bearings 

Bronze bushings, 3 A, 3/8 ID 


Shaft 

440c SS 


Link components 

7075 A1 


Mechanism weight * 

(without actuation) 

14.75 lb 

3D Mechanism weight 

(with actuation, includes a 20% factor) 

32.5 lb 

Working Volume* 

5.5 inch radius cylinder Height is 
selectable 

5.5 inch cylinder 

Stowed Volume* 

3 high x 7 radius cylinder (inches) 

7x3i inch cylinder 


lout actuators 


Analysis and design of the structural and actuation components of the 3D solar-array tracking 
mechanism are based initially on a forward dynamic load analysis tool that incorporates the inverse 
kinematics, load specifications and kinetostatic analysis to give all loads within the mechanism as a 
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function of position within its workspace. Figure 4.7. 1-5 represents some of the bearing loads and 
motor torques required in the worst-case loading scenario. The maximum of these loads was selected 
as requirements for actuator and bearing design. 


motor torques 


joint forces 
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0.4 - / \ / \ / \ 

0.2 - \ / \ / \ 
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Figure 4.7.1-5: Motor and Bearing Loads 

(For worst case loading condition) 

4.7.2 Summary of 3-5R Spatial Parallel Architecture and Kinematics 

(Pitch- Yaw-Plunge Mechanism) 

A brief discussion of the parallel architecture wrist is provided followed by two a description of two 
potential applications. 

Kinematically, the mechanism consists of 3 5-R (revolute) chains resulting in three degrees of freedom 
at the output (distal) plate. These dof allow pitch and yaw of the distal plate and translation along the 
distal plate normal axis. The 3 dof mechanism requires three inputs. This device is an in-parallel 
mechanism; there are multiple load-bearing paths connecting the output plate to the base plate (in- 
parallel). Due to this parallel mechanism topology, the device has many advantages. The advantages 
associated with this 3 dof joint include: 

• High payload-to-weight ratio 

• Good dynamic performance 

• Ground-relative actuation (all input actuators located on the base plate) 

• All revolute-joint design 

However, historically, parallel mechanisms have found a lesser degree of implementation due to a few 
limiting factors: range of motion, lack of closed-fonn kinematic solutions and lack of topologies and 
designs that match the application. In the proposed mechanism and its pointing applications, these 
limitations have been overcome as follows: 

• Closed-form forward and inverse kinematics exist 

• Large range of orientation motion (pointing) 

Finally, in addition to the positive improvements in areas of performance, there are a few key 
improvements of this approach over conventional (serial) pointing mechanisms; 

• Greater-than-hemispherical singularity-free workspace 
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• Control algorithms are not constrained by singular regions 

• Improved dexterity based on common measures of robotic dexterity. 

Figure 4.7.2- 1 below shows one prototype version of the wrist that demonstrates real-time inverse- 
kinematic control. 



Figure 4.7.2-1: Prototype 3 DOF Mechanism 

Our group has developed this device for many applications, including robotic wrists on industrial and 
EOD manipulators, high-precision manipulation devices, miniature precision control of mirrors, 
compliant, monolithic designs, and as motion simulation platforms. See Figure 4. 7. 2-2 below. 



Figure 4.7.2-2: Robotic End-Effector 

(Compliant design for pointing mirrors, sub lOnm precision manipulator) 


Two current applications are included below and in the attached document. One briefly described the 
design of the joint as a solar-array tracking device as part of the CEV, the second a gimbal mount for a 
rocket thruster (below). 


A collaborative effort between the ISPT group at MSFC, Orion Propulsion and Tennessee Tech is 
demonstrating the design of a gimbaled thruster that integrates the pointing mechanism with the rocket. 
The mount is based on the 3 dof parallel-architecture joint described above and provides the large 
range of singularity-free motion as indicated. The primary design efforts are integrating the thruster 
body as part of the distal plate and incorporating the flow-passages into the parallel-branches of the 
mechanism. Figure 4. 7. 2-3 shows a current design. 
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Figure 4.7.2-3: 3 DOF mechanism design for gimballing a rocket thruster 

Figure 4. 7. 2-4 below shows one fuel-passage incorporated into the leg. This project is also 
investigating the ability of providing real-time thrust measurements through the force-control mapping 
that is provided in the manipulator control algorithm. 



Figure 4.7.2-4: 3 DOF Mechanism with Fuel Passage Integrated into Legs 
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Figure 4. 7. 2-5 below represents the CAD models, prototype and assembly view of the Gimbal Joint 
Thruster Mechanism: 




Figure 4.7.2-5: 3 DOF Mechanism as CAD Models, Assembly 

(in hardware and prototype system) 
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5.0 DEVELOP TEST FACILITIES 

(Necessary To Replicate MXER Tether Capture Dynamics) 

This section summarizes the development of a physical capture test and the construction of a facility 
suitable for such a test. The test is used to examine the capture characteristics of a variety of 
mechanisms through the use of scale modeling. In order to simulate the dynamics of the capture 
sequence a scale model payload is propelled upward toward a scale model of the capture mechanism 
where the two will dock in a gravitationally accelerated environment. 

5.1 General Dynamics of Test and Simulated Motion of Tether 

The dynamics of the test simulate the dynamics expected of a MXER system. Due to the rotating 
motion of the tether system the tether tip traces a cycloidal motion profile thru space. In following this 
motion profile the tether tip experiences a nearly constant acceleration. When the rotational rate and 
orbital speed of the MXER system is set to dock with a payload the tether tip will briefly hold zero 
velocity relative to the payload. The extreme length of the MXER system means that only the portion 
of the cycloid which plunges in, stops, and then plunges out in the opposite direction is in the vicinity 
of the payload during capture. This means that, from the perspective of the payload, the tether tip and 
the capture mechanism that is mounted on the tether tip approach the payload in a nearly straight line, 
stop, and retreat in the direction from which they came all at a constant acceleration rate. The tests 
performed took advantage gravity to replicate the constant acceleration rate and singular direction of 
the tether motion. Taking the tether tip as the inertial reference frame the test propels a simulated 
payload upwards toward a stationary tether tip and allows the constant acceleration of gravity to pull 
the simulated payload back down. In this way the test system has a constant acceleration and the in- 
stop-out motion profile of the MXER system in a more convenient experimentation environment. 

5.2 Description of Test Facilities 

The test facility used for the capture experiments was found to be adequate and was modified to meet 
the needs of the experiment. A former racquetball court in the west side of the TTU stadium was 
chosen as the most appropriate available space for the capture test due to the size of the room. This 
space is a concrete room approximately 8m high and 6m wide. Initially the room was not equipped 
with any hardware mounting points on the walls or ceiling. In order to run the test capture 
mechanisms would need to be suspended from the ceiling so a scaffold was constructed and a 453 kg 
winch was mounted to the ceiling. Cable management hardware was also mounted to the ceiling to 
facilitate the test. Graduation lines were also painted on the walls of the room to provide some sense 
of scale in the video footage. Finally additional lighting was brought in to aid the video resolution and 
the scaffold was temporarily fixed to the wall to provide an additional measure of safety as well as to 
steady the video footage. Figure 5.2-1 below shows the resulting test facility. 
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Figure 5.2-1: Overview of Capture Mechanism Test Facility 

5.3 Fall- Arrest Pad and Safety Measures 

Safety was a primary concern during testing and a number of measures were taken to ensure that all 
those entering the test facility would leave safely. Each active component of the firing system was 
equipped with a safety off. The launcher trigger was equipped with a large red kill switch that 
prevented the system from energizing accidentally. A safety off came standard with the air 
compressor used. Every night the air compressor was decompressed to prevent an accidental discharge 
of the mechanism. A safety zone was established around the launch area and was not to be occupied 
during a test. In the event of a failed capture a “fall arrest pad” was placed under the mechanism. This 
pad was constructed of three 16cm thick foam pads which would absorb the impact of a failed capture 
preventing injury and preserving the life of the payloads. Finally a barrier ramp was constructed 
around the landing area. This barrier ramp was constructed using lxl boards and plastic netting. The 
intent of this barrier or “chicken coop” was to direct the payloads back onto the splashdown pad in the 
event that one payload took an erratic bounce. However one payload that did escape the fall arrest pad 
demolished the chicken coop and a more robust system is advised in any future testing. 

5.4 Payload Launcher Building and Development 

In order to produce the desired motion for the simulated payload a launcher device was constructed. 

At the time the launcher was constructed a target weight for a simulated payload had been identified 
how-ever a test site had not. A design goal was established which required the launcher mechanism to 
propel a weight of 4.5 kg a distance of 9 meters into the air. Because of the speed and length of travel 
needed to accomplish this goal compressed air was chosen as a power source. The necessary air 
pressure and flow rates were calculated and two 76cm long, 5cm diameter air cylinders were chosen to 
provide power to the launcher. A frame was then constructed to support two 213 cm long slide rails. 
The slide rails would mount parallel to one another so that a small platform can be bridged between 
them. Steel ropes were attached to both sides of this platform and then connected to the air cylinders 
via a block and tackle arrangement. In this way a simulated payload can be placed on the platform and 
when the system is energized with compressed air so that the air cylinders retract the platform and 
payload they are thrust upward at a high rate of speed. The entire system is actuated remotely using 
two solenoid actuated air valves. Later in the testing cycle the entire test would be fully automated 
using a LAB VIEW virtual instrument program. However details on this program are left for later in 
this report. The resulting payload launch mechanism is shown from several angles in Figure 5.4-1 
below. 
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Figure 5.4-1: Payload Launching Mechanism as Built 
5.5 Payload Fabrication and Development 

Once the launcher was constructed an array of devices were constructed to serve as simulated payloads 
for the test. Two different payload design requirements were developed. The first set of payload 
design requirements would call for a low cost, low construction time but durable payload which would 
be used for the bulk of the testing. The second payload design requirement calls for a gyroscopically 
stabilized payload which was equally durable but would be used less frequently. It was initially 
hypothesized that any design would need a considerable number of “practice shots” before the 
launcher and capture mechanism were dialed in to their preferred test settings. Therefore a cheap, 
quickly made payload that could be used for a number of “practice shots” before being discarded was 
attractive. However, at this time the test site had not yet been decided upon. The possibility that the 
payload would need to be launched as much as 9 meters into the air was still quite real. During initial 
testing of the launcher mechanism it was noted that test payloads were considerably off-pitch from 
their initial orientation on the launcher. Such a variation in the payload pitch at the apex of the 
payloads path could cause test results to be a poor representation of the real system. Therefore a 
gyroscopically stabilized payload was thought to be an attractive but expensive option. In the end an 
approach of several cheap payloads and one gyroscopically stabilized one was thought to be the best 
option available. 
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The low cost payload design was constructed as follows. Design requirements were established to be a 
nominal mass of 4.5 kg, cylindrical in shape, ~15 cm in diameter and ~3 1 cm tall. The simulated 
payload was also to be fitted with a mounting plate on the top surface which could be quickly adapted 
to a variety of booms, cameras and sensor equipment. A design was developed which made use of the 
appropriate size of schedule 40 pvc pipe, two 15cm by 4.7mm thick aluminum disks, three sections of 
8-32 threaded rod and a large cylindrical piece of brass left over from a previous project, (see Figure 
5.5-1) Sections of the PVC pipe were cut to the 15cm length and then turned on a lathe so that both 
sides of the pipe had an internal “shelf’ on the inside. The aluminum disks were then drilled for 3 
holes 120 degrees apart near the edge of the disk. Threaded rod would be run thru these holes so that 
the two aluminum disks would be clamped against the shelf in the pipe. In this way the pipe and two 
disks made a barrel shape. Inside this barrel the brass cylinder was affixed so that it would not move 
inside the payload during testing. The final payload design was of an appropriate weight, very durable 
and adaptable to many different hardware systems. Most importantly it could be replicated quickly, 
easily and cheaply since multiple collisions with the ground eventually destroyed even the most 
durable payloads. 



Figure 5.5-1: Base Payload 

Gyroscopically stabilized payloads were constructed using the following slightly more complicated 
procedure. First a section of schedule 40 PVC pipe was cut to length. This pipe was not turned to 
create a shelf but instead holes were drilled in the sides to allow switches and power recharge ports to 
be mounted. Two end caps were then turned from a 9.5mm thick aluminum plate. These end caps 
were cut so that the outer diameter matched the outer diameter of the PVC pipe but a concentric 
emboss snugly fitted into the open portion of the pipe (see Figure 5.5-2). In this way when the end 
caps are placed onto either end of the pipe the edges of the pipe are protected by the aluminum end 
cap. This design proved to be much more durable as the edge of pipe proved to be an area which was 
vulnerable to collision. The top end cap was then drilled and tapped with three 8-32 holes spread 120 
degrees from one another. The bottom end cap was drilled with three corresponding thru holes as well 
as three material relieves on the bottom side of the end cap. In this way when the nuts are fitted to the 
bottom end cap they will sink in to the materiel relief so that the simulated payload will sit on the flat 
surface of the end cap rather than resting on the nuts. The bottom end cap was also drilled so that it 
accepted a small bearing in the center of the end cap. A circular center plate was then constructed 
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from 9.5mm aluminum plate so that it slid snugly into the pipe. This center plate would also be drilled 
to accept a small bearing in the center and the three corresponding thru holes at 120 degrees from one 
another. A shaft and small commercial flywheel would be fitted to the bearings between the center 
plate and bottom end cap so that the flywheel can freely rotate between the two. The center plate 
would then be modified to accept a small 12V, 1000 rpm DC motor which engages the shaft and in 
turn spins the flywheel. It is this flywheel which will provide the gyroscopic stabilization. Since the 
flywheel is a commercially sold product it is balanced for high speeds and wobbled very little in 
operation. An 18V cordless drill battery was then fitted to the top plate to provide a compact, powerful 
and durable power supply. Finally the entire payload was assembled using three 8-32 threaded rods. 
The rods were installed so that they first screwed into the top plate. PVC pipe was then fitted over the 
rods and top plate. The center plate was then slid over the three rods and held in place by three nuts. 
These nuts were specially placed on the rods so that when the bottom plate was fitted over the rods the 
shaft holding the flywheel would be snugly fitted between the top and bottom plates so that the shaft 
was free to spin but could not jump free of the bearings. Once wired and tested this payload proved to 
be a very durable and stable test platform. 



Figure 5.5-2: Payload with Gryroscope Flywheel 
5.6 Data Acquisition Methods 

Test data for the capture experiment was recorded using a number of video cameras and 
accelerometers. Video cameras were mounted on the floor of the test area as well as at the top of the 
scaffolding to provide two separate high resolution viewpoints of the capture. A third and smaller 
camera was fitted to the tether cable at the point where the four guy wires come together. This “tether 
cam” was pointed downward at the square formed by the modified iris to provide a bird’s eye view of 
the capture. A fourth camera was mounted to the payload itself. This camera was mounted to the top 
plate of the payload and pointed upward towards the boom. In this configuration the “payload cam” 
provided an up-close view of the capture itself from the perspective of the payload. Accelerometers 
were fitted to two comer blocks. The two corner blocks were aligned caddy-corner to one another and 
each block was fitted with two one-dimensional accelerometers. Accelerometers were fitted to the flat 
faces of the aluminum angle so that they faced 90 degrees to one another in the plane of the square. In 
this way the sensors could discern the full motion of the collision between iris and payload. One 
accelerometer was also fitted inside the payload. This accelerometer was aligned vertically to measure 
the amount of shock that occurred in the payload due to capture. 
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5.7 Design and Fabrication of the Prototype Capture Mechanism 

This section describes the design and fabrication of the scale model capture concept developed by the 
TTU MXER team. The model is referred to as the modified iris and consists of a large open square 
shape which is suspended by four guy wires. When the mechanism is triggered gravity is used to draw 
the mechanism in to progressively smaller squares around a boom extended from a mock payload. As 
the squares become smaller special plates on the mechanism engage the payload boom and trap it thus 
affecting a capture. This mechanism was designed and fabricated on the TTU campus by the TTU 
MXER team. 

5.7.1 General Description of the Quad-Trap Operation 

The Quad-trap prototype is a working scale model of the modified iris capture concept. This capture 
concept uses four special booms to create a large open square at the end of the tether. During capture 
an appendage from the payload moves into the open space created by this square. The modified iris is 
then triggered and the centripetal acceleration of the tether causes the four booms to swing toward one 
another causing the square opening to become progressively smaller. As the booms draw closer to one 
another the payload appendage contacts one of the booms and is directed into one of the corner plates. 
These corner plates are mounted to special bearing blocks at the corners of the square. The bearing 
blocks slide along the booms and hold the booms together in their square configuration. As the 
payload appendage comes into contact with the comer plates it is directed into a special groove in the 
plate. The corner plates are shaped so that as the modified iris closes on a payload appendage the 
appendage will settle into this groove. When the modified iris is fully closed the corner plates 
interlock and the payload appendage is trapped inside the groove formed by the corner plates. 

The modified iris prototype works in a similar manner to the full scale modified iris described above. 
This scale model uses 4, lm long and 2.5cm diameter aluminum tubes called spreader bars to create an 
open square. These spreader bars are connected to one another via special comer blocks. Each corner 
block is equipped with two larger tubes which are arranged 90 degrees from one another and askew 
using a short section of aluminum angle. The tubes were fitted with roller bearings so that the spreader 
bars may slide inside each tube. In this way the modified iris can form smaller and smaller squares by 
sliding the spreader bars in and out thm the comer blocks. A CAD model of the Quad trap is shown in 
Figure 5.7. 1-1 below. 



Figure 5.7.1-1: CAD Model of the Quad Trap Mechanism 


Each spreader bar was equipped with a special fitting on each end. This fitting was clamped to the end 
of the spreader bar and has a small detent on the top side. The detent is shaped so that the plunger of a 
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solenoid, which is bolted to back side of the corner block, fits into the detent and holds the modified 
iris in the fully open configuration. The details of the actuation system is shown in Figure 5.7. 1-2 
below. 



Figure 5.7.1-2: Quad Trap Actuation System 

In the fully open configuration the modified iris creates an open square which is lm wide. The 
modified iris was then suspended from the ceiling via four cables which were connected to the comer 
blocks. Because these cables were connected to one another above the center of the square the 
modified iris mechanism has a natural tendency to swing closed (see model in Figure 5.7. 1-3 below). 



Figure 5.7.1-3: CAD Model of the Quad Trap Mechanism Prototype Attached to the Tether Tip 


In fact when the solenoids are activated so that the comer blocks are released from the fitting in the 
spreader bar the modified iris swings closed. However this swing was not fast enough to represent the 
behavior expected of the actual system. Therefore springs were hidden inside the spreader bars and 
cables were routed from the springs to each corner block so that the mechanism will snap shut faster 
than if it were driven only by gravity. Figures 5.7. 1-4 and 5 below show some of the details of this 
system. 
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Figure 5.7.1-4: Internal Springs to Assist Closing in Quad Trap Prototype 



Figure 5.7.1-5: Cable Management System in Quad Trap Prototype 


Finally each comer block was fitted with a specially shaped plate. This plate will serve to direct an 
appendage from the simulated payload away from the spreader bars and into a specially created notch 
in the center of the plate. The notch serves to trap the payload appendage so that it cannot escape and 
a capture occurs. These plates are shown in the CAD models of Figures 5.7. 1-1, and 2 above. 

5.7.2 Detail Drawings of Scale Model 

The detailed drawings of the scale model modified iris constructed by the TTU MXER team are 
provided in Appendix C. 

5.7.3 Autonomous Capture Functions: The Automated Capture Control System 

The automated capture control system (ACCS) consists of three primary components: a central control 
unit, a proximity sensory device, and the capture mechanism. The central control unit is comprised of 
a microcontroller unit (MCU), the MCU power supply, the capture mechanism power supply, and the 
hardware interface. The prototype configuration utilizes a Motorola FIC12 for the MCU, power 
provided from a single source (Topward 6302D) and an interface for the anning/disanning of the 
device. The proximity sensor is a Sharp GP2Y0A02YK Infrared Ranger (IR) with a range of 20cm 
(8in)to 150cm (60in). 
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The IR sensor is mounted on the PLS launch pad and positioned such that it measures the distance 
between the aft end of the payload and the launch pad. In this case the location of the payload is 
measured relative to the launch pad, but since the test setup allows for the launcher to be in a kn own 
location relative to the capture mechanism, this is equivalent to measuring the range of the payload 
from the capture mechanism. After the ACCS is armed, the MCU continuously monitors the range 
information returned by the IR sensor. This is the only external sensory device on the controller as it 
has no knowledge of the firing of the PLS aside from the range information returned by the IR sensor. 
As the payload reaches a preset distance away from the PLS launch pad, the MCU triggers the 
solenoids on the capture mechanism to release allowing the mechanism to close. After the successful 
capture has occurred, the ACCS is manually reset to await the next launch. Figure 5.7.3- 1 below 
shows this system in schematic fonn. 




Figure 5.7.3-1: Schematic Representation of the Automated Capture Control System 

Procedure for operating the ACCS: 

1 . Reset capture mechanism to open position, and payload launch mechanism to base launch 
position 

2. Place payload on the launch pad covering the IR ranger 

3. Turn on the power to the MCU for the automated capture circuit 

4. Perform a reset of the MCU by pressing the reset button 

5. Arm the solenoids by flipping the master switch to the upward position 

6. Fire the launch mechanism 

7. After capture, disable the automated capture electronics by returning both the solenoid arming 
switch and controller power switch to their off positions 

5.8 Perform Testing and Demonstration 

(On Proof-of-Concept Prototypes of Capture System Concepts Developed in Task 2.4) 

This section describes the physical testing done on the modified iris (TTU), dart (LM) and pat-trap 
(LM) mechanisms. Each mechanism was tested physically in the TTU capture lab in order to prove 
out each mechanism and ensure that no aspects of the design had been neglected. The procedure for 
conducting the tests is summarized and a brief description of the tests done is provided. An account of 
the test results is also added into this section. 
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5.8.1 Testing Procedure 

The procedure for running capture tests for the modified iris is as follows. Using the ceiling mounted 
winch the capture mechanism is lowered to a workable height. Two or more people then reload the 
capture mechanism by retracting the spreader bars to the fully open position and positioning the 
solenoid plungers into their detent slots locking the mechanism open. The payload is then recovered 
from the mechanism and reset on the launcher platform. Video recordings and sensors are then reset to 
their ready positions before the winch is used to lift the capture mechanism to its test height. Safety 
mechanisms are then deactivated before the video cameras are set to record and the launch program is 
started. The launch program uses a data acquisition board and a LAB VIEW program to fire the 
launcher, delay for a set period of time and then trigger the capture mechanism. Delay time is varied 
for each test to represent the positioning tolerance inherent in the MXER capture process. During the 
entire launch and capture process the LAB VIEW program is recording data from all five 
accelerometers for a total of 15000 samples per test at a sample rate of 500 Hz. Captures are recorded 
as successful or unsuccessful and the data recorded is stored for later analysis. Safety catches are then 
restored and the test process is repeated. 

The test procedure for the dart mechanism is as follows. A single 14-20 bolt is threaded into the 
provided hole in the rear of the dart mechanism to connect it to the simulated payload. The simulated 
payload is used with the dart mechanism because the narrow base of the device makes it somewhat 
unstable during launch. The wider base of the simulated payload helps to stabilize the device in both 
launch and flight. Use of the simulated payload also reflects a more realistic test scenario. The cables 
of the net device are checked to ensure that they are all equally tight and no individual strand has 
worked its way loose. This net device is then suspended from the winch cable and raised to test 
height. A safety procedure check is conducted before the air compressor is charged and the launcher is 
aimed for the center of the net. Air pressure is set to its highest level for this test to ensure that the dart 
will collide with the net with adequate velocity for penetration. The launcher is then fired and the 
procedure is repeated. 

The pat trap mechanism is tested as follows. Four eye bolts are used to suspend the mechanism from 
the winch cable. Drivers and motors are then connected using the diagram provided with the 
mechanism. A computer terminal is then connected to the control port of the driver board and power is 
applied to the driver as per the driver manual. Great care is taken to ensure that the current thru the 
driver is properly limited. The executable file sent with the mechanism is then opened. That 
executable file is used to slowly open the trap and hold it in the open position. Using careful hand to 
eye timing the launcher is manually triggered and the trap closed separately using the executable file 
provided with the mechanism. Safety catches are then reset and the procedure starts over again. 

5.8.2 Brief Summary of Testing Activity 

Tests were performed according to a test plan developed based on pragmatic and statistical design of 
experiment techniques. Outputs for the test were recorded as a capture/fail to capture criterion and a 
maximum acceleration shock criterion. The test plan and resulting testing processes are summarized 
on Table 1 1 below. Details of these tests are contained in the following subsections. 
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Table 11: Summary of Capture Mechanism Testing Activity 


Test 

Series 

Description 

Dates 

Summary of testing 

PLM 1 

Payload Launch mechanism 
testing 

7/04- 

10/04 

Used to perfect PLM for use on capture 
mechanism testing: Approximately 100+ 
tests performed at varying payloads 

QT1 

Initial tests of the Quad-trap 
mechanism 

2/05- 

4/05 

Tests demonstrated robust operation of 
the Quad-trap. Approximately 50 tests 
performed with >75% capture rate 

QT2 

Extended tests of the Quad-trap 
mechanism, dynamic measurement 

5/05- 

7/05 

Tests to collect key acceleration data 
during capture event on the Quad-trap. 
Data included 8+ channels of acceleration 
+ 2 channels of video. 10+ tests 
performed with data captured 

Netl 

Basic testing of the net mechanism 

6/05- 

9/05 

Four test performed before failure of 
capture mechanism. Capture was 
successful, significant dynamic transients 
occurred and resulted in mechanism 
failure.* 

BT1 

Basic testing of the Pat-trap 
mechanism 

6/05- 

9/05 

Pat-Trap mechanism tested for static 
operation (basic trap closure). Trap 
closure achieved 1 time before failure of 
drive electronics.* 

QT3 

Extended dynamic testing of the 
Quad-Trap mechanism 

10/05 - 
2/06 

Second phased of tests to collect key 
acceleration data during capture event on 
the Quad-trap. 15+ tests performed with 
data captured. Quad-trap mechanism 
required repair of pulley units prior to 
these tests 

QT4 

Autonomous operation of Quad- 
trap mechanism 

2/06- 

4/06 

Implementation of real-time ranging and 
autonomous operation of the Quad-trap 


Tests on Netl and BT1 did not receive equiva 


ent implementation period to Quad-trap mechanism by 


TTU team, and the results should be correspondingly evaluated. 


5.8.3 Test Results for the Quad Trap 

This section will briefly describe the test cases that were considered for experimental testing and 
dynamic simulation along with the data collected. First, a brief description of how the quadtrap 
prototype (and correspondingly the ADAMS model) were instrumented for testing is provided. 
Instrumentation consists of six accelerometers in three-axis triplets attached to the quadtrap prototype 
(see Figure 5. 8. 3-1) and four cameras; two cameras inertially fixed with vertical and horizontal view 
points, one camera attached to the payload and one camera attached to the quadtrap. The 
accelerometers are located at the center of mass of opposing corner nodes (see Figure 5. 8. 3-1) and are 
aligned such that the z axis corresponds to the vertical while the x and v axes lie in the plane of the trap 
and are aligned with the spreader bars. 
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Figure 5.8.3-1: Schematic of the Quad Trap for Instrumentation 

Two tests will be considered to demonstrate the capture process; 

Case 1: capture occurs in quadrant 1 (Figure 5. 8. 3-1) 

Case 2: capture occurs in quadrant 3 (Figure 5. 8. 3-1). 

Results from these tests are shown in the following section. 

5.8.3. 1 Results from Case 1 

Case 1 consists of a capture of the payload in quadrant 3. One of the inertially fixed cameras records 
the position of the payload as it enters the plane of the trap as shown in Figure 5.8.3. 1-1. Figures 
5.8.3. 1-2 and 5.8.3. 1-3 show the accelerometer data for nodes 1 and 3 respectively over a duration of 
time spanning the capture event. Figures 5.8.3. 1-4 & 5 show the individual components of 
accelerations (. x , y) for nodes 1 and 3 respectively. Each figure shows both the experimentally 
collected data along with that predicted by the ADAMS model for this capture case. 



Figure 5.8.3. 1-1: Image Record of the Case One Test 
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Figure 5.8.3. 1-2: (left): X,Y,Z Accelerations for Node 3 
Figure 5.8.3. 1-3: (right) X,Y,Z Accelerations for Node 1 



Figure 5.8.3. 1-4 (left): X (above), Y (below) Accelerations for Node 3 
Figure 5.8.3. 1-5 (right) X (above), Y (below) Accelerations for Node 1 
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5.8.3.2 Results from Case 2 


The second case consists of a payload capture in quadrant 4. The position of the payload as it enters 
the plane of the trap is shown in Figure 5. 8. 3. 2-1. Figures 5. 8. 3.2-2, 3 show the accelerometer data for 
nodes 3 and 1 respectively over a duration of time spanning the capture event. Figures 5. 8. 3. 2-4, 5 
show the individual components of accelerations for nodes 3 and 1 respectively. Each figure shows 
both the experimentally collected data along with that predicted by the ADAMS model for this capture 
case. 



Figure 5.8.3.2-1: Image Record of the Eerror in the Case 2 Test 



Figure 5.8.3.2-2 (left): X,Y,Z Accelerations for Node 3 
Figure 5.8.3.2-3 (right) X,Y,Z Accelerations for Node 1 
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Figure 5.8.3.2-4 (left): X (above), Y (below) Accelerations for Node 3 
Figure 5.8.3.2-5 (right) X (above), Y (below) Accelerations for Node 1 


5.8.3.3 Discussion of these Results 

The dynamic measurements observed in the experimental tests and those provided by model 
simulation are evaluated both quantitatively and qualitatively for the purposes of demonstrating proof 
of concept, validating the dynamic model, and providing better understanding of the capture process. 

5.8.3.3.1 Observations of the Capture Event 

The capture process that was experimentally simulated in the lab and modeled in software can be 
divided into five stages. The events that divide these stages are labeled on Figures 5.8.3. 1-2, 3. The 
first stage covers the period of time prior to initiation of the capture event. In this stage, the capture 
mechanism observes a uniform z-component of acceleration due to motion along an accelerating path 
(simulated in the lab by gravity). The second stage of the capture process starts when the capture 
event is triggered and the capture mechanism begins to close. The capture mechanism now 
momentarily moves along a non-accelerating path and the z-component of acceleration goes to zero as 
the trap begins its motion along a pendular closing path. As the rate of closure increases, the z 
acceleration increases due to centripetal acceleration about the attachment point to the tether. The x 
andy axes also observe a sinusoidal component of this acceleration (to rsin( 6)). The third stage begins 
with the initial impact occurring between the payload-boom and one side of the quadtrap. From the 
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acceleration response in Figures 5.8.3. 1-2 - Figures 5.8.3. 1-5 and the photograph in Figures 5.8.3. 1-1, 
it is observed that this first impact occurs between the payload boom and quadtrap spreaderbar #3 
(Figures 5. 8. 3-1). The fourth stage starts with subsequent impact events that may occur during 
capture; impact between the boom and a second side of the quadtrap, and final closing of the quadtrap. 
From the acceleration response in Figures 5.8.3. 1-2 - Figures 5.8.3. 1-5, it is observed that the second 
impact occurs between the payload boom and quadtrap spreaderbar #2 (Figures 5. 8. 3-1). Final closing 
of the trap is not as distinct in Fig. 8. The final stage covers the period of time after capture. During 
this time, transient vibrations die out and the capture mechanism and captured payload observe steady 
acceleration in the z direction. 

5.8.3.3.2 Comparisons with the dynamic model 

Based on this quantitative comparison of results over the series of events that define the capture 
process the model simulation is seen to correspond well with the experimental data collected. The 
order and nature of dynamic behavior over the series of capture events are in good agreement. The 
expected direction of dynamic response are consistent and the magnitudes conform (within 
approximately 10% over many portions of the capture process). Some discontinuities are also 
observed between the experimental and simulation results. First, some degree of variation in the 
timing of events is observed, such as time of the first and second impact events. This is due in part to 
the ability to estimate the error in payload position for a given experimental trial, and the inability to 
accurately represent all energy losses in the prototype system. The second variation observed is the 
degree of transient or higher-order dynamic phenomena seen in the experimental results, but not the 
model. Again, this is attributed to the imprecise nature of the prototype, in particular issues such as 
chatter, non-perfect bearings and a non-perfect inelastic collision as predicted in the ADAMS model. 
Damping was added to the dynamic model to mimic the friction introduced by the bearings and other 
such components, but an accurate measurement of the amount of damping needed has not been 
performed at this point. Third, some variation in the magnitude of z axis acceleration observed during 
the impact process. Again, this is attributed to the variable location at which the boom penetrates the 
x-y plane of the trap and the variation in vertical error of the experimental system. In some capture 
experiments, the payload is trapped while still rising relative to the trap (causing a positive z 
acceleration) and in other cases the payload is falling when capture, causing a negative z acceleration 
on the quadtrap. 

5.8.3.3.3 General observations on the capture process 

Some final observations are made on the conceptual design and the capture process. Based on the 
predicted error window and the experimental tests that were perfonned over this predicted range of 
error, the quadtrap capture mechanism proved surprisingly robust at performing the capture process. 
Capture with some degree of error (position or velocity) results in one or more impact events between 
the payload boom and the sides of the quadtrap. Fiowever, these events take place between 
components of relatively low mass (the quadtrap and payload boom) and therefore limit negative 
overall effects on the system. The acceleration forces observed on the quadtrap at times were quite 
high, at levels up to 6-7 g, but these impacts are experienced to a much lesser degree on the payload. 
Similarly, the impact forces between the payload boom and quadtrap spreader bars in some cases were 
significant, but were well within the design criteria specified for the capture mechanism. The payload 
boom in general must be the lightest possible device that could deploy in a reduced gravity 
environment, and could be damaged in the capture process. However, the payload boom has a one- 
time use and thus plastic deformation is not necessarily a critical factor. Some degree of plastic 
deformation was observed in the payload boom. It was observed over a significant number of tests 
that the system is robust to both positional error as well as timing of the release of the capture 
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mechanism. However, proper timing of the capture mechanism can reduce the dynamic impact of the 
capture process making the capture mechanism more effective. 

5.8.4 Test Results for the Harpoon and Net (LM) 

A total of four tests were conducted using this mechanism and three of those tests were successful 
captures. Tests were set up so that there was a great deal of differential velocity between the dart and 
net. In the single failed test the dart mechanism failed to penetrate the net and therefore had no 
opportunity to snag. Two of the tests snagged the net with three fingers where a third snagged the net 
with only two. The number of tests was limited to four because the fourth test broke the cable 
suspending the net and testing was discontinued in lieu of repairs. 

5.8.5 Test Results for the PatTrap (LM) 

Testing of the PatTrap mechanism consisting of static, ground-level testing of open/closing operation 
for the device. Approximately one open/close cycle was performed before a failure occurred in the 
drive electronics. The drive electronics were repaired (purchase of a new driver) for continued static 
testing. The second failure occurred almost immediately with a combination failure of the new driver 
and the mechanical mounts for one drive-train mechanism. 
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1.0 DEVELOPMENT OF THE MXER TETHER SYSTEM 


(For Dynamic Model and Algorithms) 

1.1 Overview 

The following provides discussion of an approach to tether dynamic modeling process and numerical 
analysis that derives the governing equations of a continuous elastic system representing the tether. 
Every attempt will be made to describe the scope as well as assumptions and limitations of each step of 
this algorithm. 

1.2 Equations of Motion for Elastic Body Tether 

A derivation of the equations of motion followed by a discussion of various solution procedures will 
be presented here for a continuous elastic space tether. 

1.2.1 Strain Formulation 

Consider an elastic space tether as shown in Figure 1.2. 1-1. In this model, the tether is represented as 
a one-dimensional object moving in Euclidean space. Consider the tether in both an undisplaced and 
displaced condition as shown in Figure 1.2. 1-1. 



Figure 1.2. 1-1: Tether in Displaced and Deformed Condition 



Figure 1.2. 1-2: Displacement of Tether Element 

In this figure, two metrics for parameterization are shown, s, a material coordinate and 5*, a geometric 
coordinate. The coordinate s is invariant with respect to space and time while the geometric 
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coordinate s * is uniform with respect to geometric or defonned length. Due to static and dynamic 
forces on the tether, the tether undergoes a combination of rigid body and elastic displacements. The 
displacement of a small element from ds to ds * is shown in Figure 1.2. 1-2, where u, v and w represent 
displacements in the x, y and z (not shown) directions respectively and a trailing prime indicates partial 
derivatives with respect to space (geometrical space coordinate). The geometry of the displaced 
element can be defined as shown in Figure 1 .2. 1-3. 




Figure 1.2. 1-3: Tether Element Geometry 

The length of this tether element can be defined as, 

ds* = ^(x'+u'f +(y'+Vf +(z'+w'fds 1 

where again s is the coordinate based on the tether material, while s * is the coordinate based on current 
tether length. From these equations, the total length of the tether can be determined at any time as, 

/* = j' = ( ' ^{x'+u'f + ( v’+v’) 2 + (z'+w'fds . 2 

The strain can be defined at any point along the tether as a function of tether element displacement. 

An engineering strain definition will first be implemented for the tether as, 

ds * -ds _ds* i 
ds ds 

with s the strain given as a function of space and time. The ratio of stretched to original element 
ds ^ I 

length, — — , will be defined as rj with r) = tJ(x'+u ') 2 + (v’+v’) 2 + (z'+w’) 2 . r\ gives a relation between 
' ds 

the element differentials, ds = rj ■ ds . 



1.2.2 Reference Coordinate System 

First consider an reference system based on the coordinate, r(s,t) which will represent the tether 
position with respect to an inertial reference frame (selected as an earth-centered, non rotating frame, 
(Figure 1.2. 2-1), defined as; 

r (s, t) = fx + u)i , (y + v)j, (z + w)k,\ 
r’ {s,t)= | x'+u')i , (y’+v’)y , ( z'+w')k , 


Eq. 4a 


or 

r x ' = (x'+u'\ r v ' = (y’+v’), r, ' = (z'+W) 4b 

With r = r (s,t) the primary parameter describing the tether with respect to space and time, the strain 
equation can be written as: 

s(s,t) = r]-l 


with 
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n = j{r x 'Y +M 2 + (c ') 2 


4c 



Figure 1.2.2-1: Tether System 


1.2.3 Tether Geometry 

This section will describe the process of making transformations between the tether frame and the base 
(fixed) frame or coordinate system 1 . Consider an inertial base frame, {B} (x,y,z), and a tether-fixed 
frame, {T} (u,v,w) such that u(s,t) is the tether tangential direction at s, v(s,t) is the tether radius of 
curvature direction at s and w(s,t) defined to complete the right hand frame (fig. 5). The tether 
tangential is given in r’ (called x) while the tether radius of curvature is given as (r” = x’)- The 
magnitudes of x and = x’ are defined as // and k respectively 

( rj = J Xx + X y + Xz ■> K = Jx'x + X'y + X'z )• Then, for the general case of tether orientation in 

space, the tether frame orientation with respect to the base frame is defined by the rotation R T b (see 
Figure 1.2. 3-1 below) This rotation matrix (defined to be a member of valid homogenous rotation 
matrices) is given as, 


Certain parameters will be easier to describe first in the tether frame, such as internal tether forces/moments and tether mass moment of inertia, and 
then converted to the base frame (all motion must be described with respect to an inertial or Newtonian frame for later developments). 
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Figure 1.2.3-1: Tether Frame (u,v,w) and Base Frame (x,y,z) 


The rotation from fixed to tether frame can equally be defined as the Euler (&,-</>, y; w, v, it) rotation 
set, 


cdctj) 

-sOcy -cOs(j)sy 

sQsy -cOs(j)cy 

sQscj) 

cGcy -sQscj)sy 

-cOsy -sOscjxzy 

S(f) 

ccjtsy 

ctjtcy 


The derivation is readily seen based on the fact that projections of the tether tangent onto the fixed 
frame are (r’ x , r’ v , r’ z ) and the projections of the tether curvature direction onto the fixed frame are 

(X x, X y’ X z)- 

For this geometry, a number of special cases may exist. For example, the case of planar motion, or 
motion in which the tether tangent and curvature vectors lie in the orbital plane (with x,y defining the 
orbital plane) is considered. In this case, angles (j) and y in equation 5b are both zero, allowing the 
rotation matrix to become, 


R 


T 

B 





0 

0 

1 


6 


An example of use of this homogenous transformation is provided. Consider the reaction forces, P at a 
point along the tether consisting of an axial (tangential) force, N in the u direction, an in-plane shear 
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force acting in the direction of curvature (v) as Q, and an out of plane force O acting in the w direction. 
This reaction set can be resolved into the base coordinate system as, 



1 

A, 

1 

. v. , J-S 

■ -sy 

o 

i 

N 


r' / v* / 

Ny -Q V 
/ rj ^ / 77 

P B = R B P T = 

r / r / 

y V o 

A A 

0 0 1 

0 

0 

> = < 

Y / y ’ / 

N V +Q y 
/r/ ^ /rj 

0 


Figure 1.2. 3-2 shows the tether with base and tether frames attached. For purposes of this work, the 
inertial frame {B} can be given as any convenient “near inertial” frame, such as an ecliptical 
coordinate system or equatorial coordinate system. 



Figure 1.2.3-2: Tether System with Frames 
1.2.4 Equations of Motion from Newton’s Second Law 

A free-body diagram (fbd) of an element of stretched length, ijds of the tether is shown in Figure 1 .2.4- 
1 . This element undergoes a combination of distributed and concentrated forces. 



Figure 1.2.4-1: FBD of Tether Element 
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The forces acting on the element are summed and set equal to the time derivative of linear momentum, 

(pAds dr / dt) = P + ^/q s ds - P + ¥ d ijds + F, 8a 

or 

pAr = P' + F 77 8b 

which assumes that p and A are functions of space but not functions of time, r, P, F d , F d , and F| are 
functions of space and time and an over dot represents derivatives with respect to time and a trailing 
prime representing derivatives with respect to space (the material space coordinate). Note that F 
represents the resultant of the external forces given on a per length basis (F = F d + F|/(ip/.s')). These 
forces will consist of a combination of gravity, electrodynamic drag, aerodynamic drag and the like. 

In a similar manner, the sum of moments acting about one end of the tether are evaluated, 


- dr / dsx(pAdsd 2 r/dt 2 )+dIa + (oxdIo) = (M + dM / ds)-M + dr / dsxl P + 5I % ds)+- dr /, dsxFds 
2 /as ' /as > /as ' /as > 2 '® s 

9a 


or 


^-r 'x[pAr)ds 2 + dla + to x dla = M 'ds + r'ds x (p + P'c/.y)+ r' x F(.s’, t)ds 


9b 


where parameters r, P, F, M, I, a, and co remain functions of space and time and p and A are functions 
of space only. Note that in Eq. 9 that the resultant force is assumed to act at the center of the element. 
Equations 8 and 9 represent general equations of motion for the tether element. 


1.2.5 Forms of Equations of Motion 

The following sections will provide several reduced forms of these equations. The following cases 
will be presented. Case 1) will consider tangential loading only within the tether (tether carries no 
bending or shear) and will assume the tether mass to be lumped at the element center. Case 2) will 
extend case 1) and will consider tangential as well as shear and bending loads, all occurring within an 
orbital plane (each end of the element considered to have three dof). Case 3) will add the potential for 
torsion to case 2 as well as consider mass moments of inertia of the element. Case 4) will consider a 
case of spatial force and bending (six dof at each end of the element). 


1.2. 5.1 Case 1) Tether with Tangential Loads Only 

The equations 8 and 9 are first reduced first by dropping all higher order space derivative terms, and 
assuming the element such that dl is negligible. This yields, 

pAr = P' + Fr/ 10a 

0 = M' + r'xP 10b 


In case 1, the tether is assumed to carry axial loads only (behaves as a string). In such a case, the 
element reaction end load P is given as, 

P = Mi 

where N is the magnitude of the tensile force and u represents the tether tangential direction (the 
direction of r’). Equation 10b then vanishes, while Eq. 8b can be written in scalar form for the base 
coordinate system, x,y,z using the rotation operator of Eq. 5a as, 


pAr k = 


/ r 
n’a- 


V 


V ) 


+ F k r/ , k = x,y,z. 


where N, the tensile force is defined from the tether material constitutive relation as, 

N = EAs = EA(tj -!) 


11 


12 
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where E and A are functions of space, and s and i) are functions of space and time. Thus Eq. 1 1 
becomes: 


pAr k = 


EA (// - \)r 
V 


+ F k rj, k = x,y,z. 


with the function 77 given in Eq. 4c as, 


V = jr x ' 2 +r' 2 +r z ' 2 


13 

14 


Equation 13 represents three, second-order partial differential equations (pde) in space and time that 
can be solved in terms of the unknowns, r x (s,t), r y (s,t) and r z (s,t). The solution to this boundary value 
problem requires initial conditions of the state and two boundary conditions to be specified. A 
discussion of general solution approaches for this pde is discussed in a later section. 


1.2. 5.1.1 Expand Partial Differential Equations of Motion 

The partial differential equations developed in Eq. 13 are expanded to yield the following set of 
relations, 


pAr k = EA[ifr' k + 77(77 - 1 )/;"]// 2 + F k q ,k=x,y,z 


15 


with, 


n' = [« 




+ r.r 



15a 


In this expansion, p and A remain functions of s, but A ’ is assumed to be 0. Equation 15 clearly shows 
the order (second order in space and time) as well as the nonlinearities of this problem. 


1.2. 5.1.2 State-Space Partitioning 

A state-space representation of these partial differential equations of motion is constructed to result in 
a 3» first order equations of motion (where n is the number of second order equations of motion). The 
third set of equations is based on the time derivative of % . This state-space vector is defined as, 



r 


r 

x = < 

V 

> = < 

dr/ 

Zdt 


X 


dr/ 

. Zds, 


The equations of motion can then be written as; 



r 



V 

X = < 

V 

> = < 

1 

pA 

(EA(tj - l)r'/ 77 ) +F ?7 


X 


v' 


where the third set of relations is derived from the assumption that tether motion is separable in space 


dr 

and tune, or 

d 2 r 

= . The expanded relations in Eq. 15 can be written as, 


dsdt 

dtds 


with, 

pAv k = EA[ri'% k + rj k (tj - \ )x'}j 2 + F k ri , k = x, y, z 

18 


/ = U-rl'r +Z y Zy +XzX'z\l * 

18a 


/ 2 , 2 , 2 
7 = VXc +Xy +Zz 

18b 


1.2.5. 1.3 Boundary Conditions 

Two boundary conditions are required for the second-order pde tether model. In the MXER tether 
model, the tether boundaries consist of end-bodies whose dynamics are coupled with those of the 
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tether. Therefore, the boundary conditions are defined from the equations of motion of the end-bodies. 
Consider first the case in which the end-bodies are attached to the end of a fixed-length tether. A firee- 
body diagram of a general end-body is shown in the following Figure 1.2.5. 1.3-1. 



Figure 1.2.5.1.3-1: FBD of End-Body on Tether System 


ma b, = F « + P W 

BC. la 

Thus, the tether end-load is given as -Pbi at s=0, / or, 

«M a M = F h +N(0,t)n(0,t) 

BC. lb 

m b2 a b2 = ¥ b2 ~N(l,t)u(l,t) 

BC. lc 


Next, consider the case in which the tether length is varied by reeling tether either in or out of the end 
bodies. In this case, the boundary conditions must account for this state of mass flow of tether, as 
given in the following equations, (see Figure 1.2.5. 1.3-2,) 



Figure 1.2.5.1.3-2: FBD of End-Body on Tether System with Tether Reeling 


in 


or 


w a M =M0M0M0,t) + Y hl+ N(0M0j) 


«H a M = ' 


pAih 


rj 


+ F m + 


EA{rj- l)x 


(o,0 


77 


l(o,0 


m b2 a b2 ~ 


pdl 2 X 


V 


+ F 62 + 


EA(tj- \)i 


(/. 0 


V 


m 


BC. 2a 


BC. 2b 


BC. 2c 


1.2. 5.2 Case 2) Tether with Bending in the Orbital Plane 

Case 2 considers the tether with both axial and bending stiffness, all acting in the orbital plane. In 
such an example the element end reaction P becomes, 

P = Mi + Q\ 

while the end moment M’ becomes 

M' = MW. 


87 



Now, starting with Eqs. 10 (EOM with higher order differential terms dropped), 

pAr = (Mi + Q\) + F ?7 Eq. 19a 

0 = M'w + i]ux( y Nu + Q\) Eq. 19b 

describe the modified EOMs where certain elements are described in the tether frame for convenience. 
Eq. 19b is expanded to yield, 

0 = (M' + r/Q)w Eq. 20a 


i - _ M '/ 


Eq. 20b 


Assuming a simple beam theory for the small tether element, the bending moment is proportional to 
curvature and bending stiffness as, 

M(s,t)= El(s)\r"(s,t)\ . Eq. 21 

So Eq. 20b becomes, 

El\x'\ 

Q = — • 22 

V 

Eq. 22 is substituted into Eq. 19a and transformed to the base coordinate system to yield the system of 
equations, 

^ y h ( k’ Y ^ 

pAv x = EA{rj - 1)— + El— — +F x ?i 23a 

l n ) \ n n ) 

r V 3 f k’ Y \ 

pAv = EA{rt-\)— - El — ^ + Fn 23b 

' l n V n n ) 


pAVj = EA [// ’ x, +(tj- 1 )nz'j )l 2 + El[rj(K"z k + i<x' k )-^i<<X k }l ’’ + F jij , 

(j,k) = (x,y), (y,x), 24 

with 

K = ^x'l+x'),+x') 25a 

K ' = Iz’xZ'x + X'yX'y + Z'zZ'z V' - 25b 

= l-kz: +z'yz; +z:z:k+lz:z:+z:z: +...kk- 2 . 25c 

Thus it can be seen that the new system equations (24) are second order in time and fourth order in 
space. For this system, four boundary conditions are required. The boundary conditions will be 
defined as the end force and moments that satisfy the equations of motion for the end bodies. 


1.2.5.2.1 Boundary Conditions 

The boundary conditions for case 2 consider both tangential forces and moments that satisfy the 
equations of motion for the end bodies. In this case, assume that the tether is attached to the end-body 
through a pinned connection (frictionless). In this case, the additional two boundary criteria require 
that no moment can be passed through the end of the tether. Hence, the boundary conditions become, 


m M a bl = ~ 




+ F M + 


ea{j] - !)/ 


BC, 2b 


m b2 a b2 




+ F Z>2 + 


EA(tj- l)x 


BC. 2c 
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BC. 3a 


Xl ( ,,)=° BC3b 

1.2. 5.3 Case 3) Tether Element with Mass Moment of Inertia 

In all prior cases, the assumption that the tether element can be treated as a particle has been made, 
leading to an elimination of dl terms from Eqs. 8 and 9. Case 3) will modify Case 2) by including the 
mass moment of inertia tenns. The justification for the removal of rotational inertia terms in previous 
cases was made on the assumption that a given tether element would remain small. However, in 
solution, it may be desirable to reduce the overall number of tether elements. Therefore, the ability to 
account for moment of inertia effects may be desirable. Rotation of the element is defined almost 
completely from the end node motions through kinematics, an extra degree-of-freedom (dot) is 
included to account for rotation about the tether element centerline. The kinematic relations that 
define the angular rate are given in the following equations (derived from the concepts of 
instantaneous screw axis). The result is a functional relation between angular rate, to and nodal 
positions and velocities, r and v. 

A V / X V i+1 

OJ = r 3U 

K-xv,. +1 | 

H r i - r o| = V , 

H r i + i- r o| = v ; +1 31 


®N' r o| = v ^ 

Equation 30 is used to detennine the direction of the angular rate vector, while Eq. 3 1 is used to 
determine the location and magnitude of this vector (magnitude given by co, location given by r 0 ). The 
angular acceleration term, a, can be determined from a difference approximation of a=dco/dt. The 
mass moment of inertia of the element is 

I = Jr 2 dm = jr 2 pAds . 32 

With this information, Eq. 9b is reevaluated. In evaluating the moment of inertia related tenns in Eq. 
9b, the body fixed frame will be used to perfonn the initial vector manipulation, and then as a last step 
converted to base frame coordinates. Also, it is noted that as a matter of convenience the body fixed 
frame will be assumed to represent the principle frame for the tether element allowing terms in Euler’s 
equation to reduce as follows, 

f/ ] [\I -I jco co ) 

uu u \ ww yy / ^ v w 

Ro \l do \ + \{l -I \o co \ . 33 

B vv v V uu ww / u w 

I co [I —I )co CO 

When these equations are substituted back into Eq. 9b (with higher order differential terms dropped), 
the equations become, 

4.A] JG0 u 

' 7 V>V [ + j (luu ~ L H [ = - |x | 0 34 

I zz d> z 

From Eq. 34, shear terms O and Q can be solved and substituted into the translational dynamic 
equation. The additional degree of freedom, rotation about the tether axis, defined by q u will be 
included in the state space representation of the tether. Euler’s rotational equations will be employed 
to give a first order differential relation between 9 U and co (using the Euler rotation operations 


KafiK 
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discussed in section 2.1). Thus, the state space system becomes a system of (3(3)+l)*« equations that 
can be solved forward in time. 

1.2.6 Application to Variable Length Tether 

The equations of motion resulting from the case in section 2.5.1 were given in Eq. 17 as, 



r 



V 

Ed) = < 

V 

> = < 

1 

pA 

(. EA(?i - l)r 7 tj ) + F 77 


X 


v' 


For a fixed-length tether, the spatial domain for these equations is, 

0 <s<l 2 

where the length of the tether is initially defined and fixed in time. This case allows direct spatial 
discretization to be applied (for example through a finite difference or finite element approach) with 
fixed boundaries on the domains. Alternatively, a variable length tether model may be desired, for 
example if tether reeling is to be permitted at either boundary. In this case, the spatial domain of the 
tether becomes dependent on time as, 

0<s< l(t ) 3 

where, 


t=t. _ . 

lit ) = /q + J* (/j + i 2 jdt 


t = o 


with la the initial tether length and /, , U the rate of tether reel-out at each end of the tether. In order to 
maintain constant domains for future spatial discretization procedures, a new spatial variable E, is 
introduced (following the procedure and notation introduced in Kohler, [zz]) to replace 5. ^ will be 
defined as, 

£ = - 5 

/ 

which has the fixed domain of 0 -> 1 . Based on this new spatial parameter, we can define the tether 
coordinates with the function, 

r(£,f) = r(s,f). 6 

which is dependent on variables E, and t and is represented with an over tilde. Based on this new 
parameterization, the derivatives of the tether coordinates become, 


f dr 

Xk= r k 

dE, 

dr \ l dr k 

v ‘‘ r ‘ = a + ^ 


k ~ 'k 

5 2 - 


dt 2 dE,dt d% dtd$ d % 2 


with derivatives of the new spatial coordinate given as, 




ih +0 


| = -4/+4/ 2 - 

i 2 r 


7 

8 
9 

10 

11 

12 


and the stretch ratio defined as, 
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CFV (drV 


77 + {d{} + Ud ■ 13 

The remaining state variables (and stretch ratio) are defined as functions of the new spatial parameter 
as indicated by an over-tilde as, 


n = Jz x 2 +z 2 +z _- 2 


Allowing the variables from Eq. 1 to be written in terms of the new state variables, 

Z k = ?Z k 

y < = e>2?Zk 

k /“v C* 


v k= v k+GZk 


y l= ^L + 2 ^ + S, + i^ 

k dt dg k 

v = 4'v 

eg ' eg eg 

The equations of motion in terms of the new parameterization become; 


dt pA y eg ' eg eg J eg 

= 18 
dt 


with the new state vector. 


k = x, y, z 

H-vl 

Xfe()= vfel) 


a function of the new spatial variable £ and time with domains, 

0 < t < oo 

0<£<1 ‘ 


1.3 Compare with Bead & Spring models 

Consider the bead and spring model for an elastic tether in two-dimensions with axial stiffness only. 
Figure 10 shows one section of tether between two beads. 
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JO 



The tensile force s between 


the beads is given as a function of tether stiffness and bead positions as, 



Expanding, 



+ 


r y + / 



V 


/ 


3-2 


or 



with higher order terms neglected, 5 becomes, 


s 



1 


3-3 


3-4 


with k = . Thus, the bead and spring model is consistent with Hooke’s law relating Eulerian 

stress to engineering strain. 


1.4 Alternative Reference Coordinate System 

An alternative reference coordinate system for describing tether motion is considered. Numerical 
solutions of the equations of motion, for example Eq. 18, are used to predict tether behavior. An 
alternative reference system for the tether is suggested as one approach to improve or validate the 
numerical solutions. The initial system tracks each point the tether relative to a fixed position (for 
example earth center), resulting in large numerical values for position coordinates relative to the other 
parameters. The first alternative considered here will describe the tether motion relative to another 
orbiting body. In this case, the orbiting body will be an “equivalent” rigid tether model, equivalent in 
terms of end mass and tether mass. The tether state will now be composed of the position and 
orientation of the rigid tether (6 dof) plus the three position coordinates of each point along the tether 
relative to a point on the rigid tether. The total number of coordinates is given as 3*n + 6. 
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A point on the tether is referenced with vector r as, 

r = r c +s + u 4 

and in-plane coordinates, 

r = r a + s cos(9) + u x 
r y = r cy + s sin(#) + u y ' 

The strain follows from the partial derivatives of these coordinates with respect to the material 
coordinate, 

dr/ ( n \ , du r / 

Yds = cos (9+ Yes 4 

dr / . (n \ du v / 

y ds = un(e)+ >/ ds 

/ \ ds * -ds ds* . 

s{s,t) = = 1 = 77-1 4 

ds ds 


+ K 


= 1 + 2 ( 008 ( 6 * )u[ + sin (0)u' y )+u' x 2 + u' y 


Here, the deflections ux and uy are specified in the base coordinate frame and are functions of space 
and time. These deflections describe the location of any point along the tether with respect to the 
rigid-body tether. Note that these deflections account for all departure from rigid body tether motion 
(elastic deflections, dynamics, etc). For convenience, these deflections will be written in the rigid 
body frame, u and v as; 

u x c6 - s6 0 u 

<u y >= s6 c6 0 < v > 4-6 

0 0 1 w 


u B = Ru t 4-7 

with B, T referring to the base and tether frames respectively and cO, s 6 *as cos($), s i n (T9) . Similarly, 
the partial derivatives of these deflections are given as, 

u B = Ru t 4-8 

allowing the strain to be rewritten in tether coordinates as; 

s = 77 - 1 4-9 

with 

77 = V(l + w ') 2 +v ' 2 4-10 
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Similarly, the position can be rewritten in terms of deflections in the tether frame as, 

r =r +sc0+uc0-vs0 

. 4-11 

r = r + ss0 + us0 + vc0 
Expressions for the tether acceleration then become, 

r x = r cx -(s + u)0s0 -(s + u)0 2 c0 + iic0 -2u0s0 -vs0 -2v0c0 -v0c0 + v0 2 s0 

r y = r vx + (s + u)0 c 0 - (s + u)0 2 s 0 + ii s 0 + 2u 0 c 0 + v c 0 - 2v0s0 - v0s0 - v0 2 c0 

These expressions are then placed into the equations of motion for a tether with negligible bending 
stiffness (Eq. 13), 

J EA iv~ l ) r A , r 


pAr k = 


+ F k rj , k = x,y,z. 


pA\r cx - (5 + u)0 s 0 - (5 + u)0 2 c 0 + ii c 0 - 2u0 s 0 - v s 0 - 2v0 c 0 - v0 c 0 + v0 2 s 0 } = 
' EA[r) - lX(! + u')c0 - v'sff) ^ +p ^ 

v 0 J 

and 

pA \r yx + (5 + u)0 c 0 - (5 + ii)0 2 s 0 + ii s 0 + 2u0 c 0 + v c 0 - 2v0 s 0 - v0 s 0 - v0 2 c = 

EA{r/ - lX(! + u')s0 + +frj 

v V J 

for the in-plane coordinates. These equations described in the tether frame coordinates, 


4-13a 


4-13b 


pA \r cx c0 + r s 0 -{s + u)0 2 + ii - v0 - 2vffl = [ EA ^ — ^ + “ ^ + F x tjc0 + Fns0 4 

\ 0 ) 

r 

pa{- r cx s0 + r cy c0- (5 + u)0 -2u0 + v0 2 - v}= EA ^0 — L)l_ - F x rjs0 + F y rjc0 4 


4- 14b 


ii = — \ea{i]'{\ + u ' )q 2 + (77 - 1 )u " // 1 )+ (f x c0 + F y s - \r cx c0 + r cy s 0-{s + u )0 1 + -v0 - 2vS\ 

pA 

v = — \EA(rj'v'rj 2 + (77 - \)v"rj 1 )- (f x s 0 + F y c 0)q }- {- r cx s 0 + r cy c 0-{s + u)0 -2 u0 + v0 2 } 
pA 

4-15 

These equations are coupled with the rigid body equations of motion, 
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o -[m 2 + m y^u(y-o) 

(m, +m 2 +m 3 )r ( 2 j \n 2 + m yQr c lc(r - 0) , 

- ^ m 2 + m y/ j/ ^ m 2 + m y / jr./c(/ -0) ^ m 2 + j 

- (/«j + nz 2 + m 3 )r c y 2 + (m 2 + m 3 l2il8{2f-8)i(r-8)) ]’ 

_ 8 ,H /‘ ^ ^ + l c ^y _ ■- _ pj rj ^ r - + c ^y _ 0^p / 

(m l + m 2 + m 3 )2r c r c y + (m 2 + m 3 /2)r c l 6 1 s(y - d) 

+ /jm 2 r c l s(y - 6)a % + jum 3 r c j/^ s(y - 6)b % 

(m 2 + m 3 /2ir c yl(y - 2&)s (y -&) + 2r c lGe.(y - $)) 

- jum 2 rJ s(y - 0)a ^ - jum 3 r c jy^ s(y - 6)b ^ 


A = 


(m ] + in 2 + m 3 ) 
0 


4- 17a 


4- 17b 


with r c , /locating the reference position on the rigid tether and <9 the rigid tether absolute orientation in 
the orbital plane, ml, m3 the end body masses, m2 the tether mass (tether mass assumed to be 
uniform), / the tether mass moment of inertia, / the rigid body tether length, and 

[(/; c A + / c 0) 2 + {r c s A + / s 0) 2 ] 


a = 


b = 


c A + 


y 2 c8 


■M 


r„ s A + j/^ s 6 


4-17c,d 


This set of equations of motion, 4-15 and 4-16, are equivalent to those given in Eq. 13 and can be used 
to evaluate the numerical solutions imposed. An example state vector for the equations of motion from 
section 4 based on a finite difference solution approach with n nodes might look like, 

V 

7 

5 




4-18 


1.5 Tether Damping, Strain-Rate-Proportional Damping (Material) 

A model that includes strain-rate proportional damping is provided here. This damping is included in 
the constitutive relation (Eq. 1 1) as, 

N = {EAs-cs ) 5-1 

where c is the strain-rate damping coefficient and the strain rate defined as, 

^ = |( ? 7- 1 )= 7 7 = U v Xx + XyXy + XzX:\l 1 • 5-2 

The equations of motion that include this strain-rate proportional damping become, 
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pA\ = 


5-3 


EA(r/-l)+cfj 


V 


X 


+ F . 


Upon expansion, this system of equations in scalar form becomes, 

pAv k = [EA(ri'z k + )J k (>? ~ 1 )x') + c (>l>lX k + VVX' k ~ VV'X k )]>1 2 + F k , 

k = x,y,z 

with 

/ 2 ! 2 ! 2 ~ 

V ~\X x + Xy + Xz 

1)' = iXrX'r + XyXy + XzXz}l ' 

V = U, j.v + Zv/fv + Jz Xz \l X 

n = l(x'X x +X x X' x +XyXy +XyXy + X'zXz + X X z \l ~ U X X x +XyXy+X-_X : )n , Y 


5-4 


5-5a 

5-5b 

5-5c 

5-5d 


2.0 VALIDATION OF TETHER SYSTEM MODELS 

In developing any engineering model and subsequent computational solution schemes, validation 
forms a critical process. Many traditional techniques for solution validation exist including 
convergence tests, comparison with closed-form solutions, and tests with physical experiments. There 
are many situations however when some of these techniques cannot be directly applied, for example in 
comparing a solution with direct physical evidence. Such situations often arise in the aerospace and 
space technology arenas and demand therefore a greater emphasis on validation of engineering models. 
This chapter will demonstrate a specific model validation procedure that will be applied to multiple 
formulations to the tether dynamic simulation code. 

MXER Tether simulation benchmarks and validation/e valuation of specific tether model algorithms is 
perceived to be a critical element in the MXER system design and development process. This section 
will discuss activities toward model validation through a multi-stage process that includes, 1) model 
validation against simple lab tests, 2) model validation against limiting cases pennitting approximate 
closed-form solutions, 3) model validation against previous tether missions, 4) verification across 
multiple alternate formulations and 5) verification across multiple, separate parties involved in model 
formulation. The initial goal of this work would provide a benchmark solution set, future work will 
evaluate the efficiency of specific model formulations in achieving a desired precision. 

The equations that govern tether behavior are invariant to method of model derivation with respect to 
the underlying physics. However, the interpretive procedures used to describe these physics are often 
approximate or not fully understood. Further, these equations result in differential functions of space 
and time which in the most general sense do not permit exact solutions. As such, approximating 
techniques are applied to these governing equations to arrive at the desired result, prediction and 
analysis of the dynamic tether behavior. In general, these approximating techniques can be divided 
into spatial approximations and temporal approximations. The spatial approximating techniques are 
generally applied to reduce the partial derivative equations (pde’s) into ordinary differential equations 
(ode’s) or into a form that allows a specific, known solution. The temporal approximations provide 
numerical integration procedures to solve the time dependent derivative functions. Taken together, 
space and time approximating techniques create a process to solve the governing tether equations. The 
techniques for space and time in general can be treated independent of each other and allow a variety 
of combinations. Further, these techniques provide some level of variation and should find 
applications in which they are more ideally suited. Therefore, a process of examining a variety of 
these techniques will both provide a process they may indicate convergence toward a benchmark 
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solution to the inverse dynamics problem as well as indicate formulations that are most ideally suited 
to this problem. Table I provides an overview of spatial and temporal techniques that are being 
implemented by TTU: 


Table 1: State versus Time 

(Approximate techniques that could be considered) 


Space / Time 
Approximation 

4 th Order 
RK, 

6 th ORK 

Adaptive 

time-step 

routines 

Adams- 

Bashforth 

Moulton 

Explicit/Implic 
it Adams 
Moulton 

Implicit 

Gear 

method 


(RK78) 




Finite difference solution 
(over multiple coordinate 
systems) 

X 

X 

X 


X 

Assumed modes based on 
Ritz formulation 

X 





Finite element formulation 






Discrete finite element 
formulation 

X 






2.1 Procedure 

The validation process will proceed through a sequence of test cases prepared for comparison either 
with some known or predicted result, or for comparison among multiple solution techniques. As such, 
the cases will be classified as either type A (known or predicted results and over alternate numerical 
strategies) or type B (comparison with independently developed models). Table II presents an 
overview of a series of cases that will be involved in the validation process and provides an overview 
of the current state of the validation process. Discussion of the results from Table II follow. 
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Table 2: Validation Cases versus Model Formulation 


Validation 
Cases -> 
Model 
formulatio 
n | 

Mods. 

Case 1) 
(type A) 
Compare 
with the 
Wave 
Equation 

Case 2) (type A) 
Spatial/tempora 
1 convergence 
models 

Case 3) 
(type A) 
Compare 
RMA over 
time 

integratio 
jn routines 

Case 4) (type 
A) 

Compare RMA 
over spatial 
approximation 

s 

Case 5) 
(type B) 
Symmetric 
, Non- 
symmetric 
tether 
rotating in 
free-space 

Case 

6) 

(type 

B) 

RMA 

Finite 

difference 

solution 

(inertial- 

coordinates) 

RK4 

Sect. 2.2.1 

Sect. 2.2.2 

Sect. 2.2.3 



Sect. 

2.3 


RK6 






Sect. 

2.3 


ABM 



Sect. 2.2.3 



Sect. 

2.3 


3-Pt. 

FD 




Sect 2.2.4 


Sect. 

2.3 


5-Pt. 

FD 




Sect 2.2.4 


Sect. 

2.3 


RK78 





Sect. 2.3 

Sect. 

2.3 

Finite 

difference 

solution 

(rigid-body- 

relative 

coordinates) 






Sect. 2.3 

Sect. 

2.3 

Finite 

element 

formulation 

(inertial 

coordinates) 








Finite 
difference 
formulations 
from TSIM 






Sect. 2.3 

Sect. 

2.3 

Assumed 
modes based 
on Ritz 
formulation 









CFS - closed-form solution 

RMA - representative MXER application 
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2.2 Results of Type A comparisons 

This section contains a description of the results of specific tests performed as indicated in Table II. 
Table III below contains the specific data used to initiate the representative MXER application 
containing two design nodes. The performance over a single orbit is presented in the subsequent 
graphs (approximately 1 100 seconds). 


Table 3: Parameters for Representative MXER Application 


Constants 


Tether Facility Orbit 


F 

3.986E+14 

Perigee (m) 

6,785,136.00 

Number of Design nodes 

2 

Apogee (m) 

14,823,136.00 

Design node 1 mass 

10000 kg 

Time for one orbit (s) 

11,176.25 

(ballast mass) 




Design node 2 mass 

500 kg 

SMA 

10,804,136.00 

(capture mass) 




Tether Length (m) 

90000 

e 

0.371987172 

Cross Sect Area (m 2 ) 

1.9635E-05 

Perigee Velocity (m/s) 

8977.690593 

P (kg/m 3 ) 

970 

Apogee Velocity (m/s) 

4109.444293 



co of tether fac. (rad/s) 

-0.013881192 


2.2.1 Case 1 

Case 1 consists of validation of the finite difference code model (4 th order Runge-Kutta integration and 
3 point finite difference routines) against an exact solution for the vibration of a simple string. The 
string is massless, I s fixed at both ends, has parameters chosen to provide for the first natural 
frequency, and is given an initial starting condition that matches its first mode shape. 

Results: 

The model was compared using three scenarios. The first case uses 20 nodes and a time step of 0. 1 
seconds. The second case uses 20 nodes and a time step of 0.05 seconds. The last case uses 100 nodes 
and a time step of 0.05 seconds. These results are shown the following 6 Figures. Figures 2.2. 1-1 and 

2.2. 1- 2 are representative of n = 20 and At=0 . I for one single node (middle node) along the tether. 
Figure 2.2. 1-2 shows the error over time between the exact solution and the model for the same node. 
Similarly, Figures 2.2. 1-3 and 2.2. 1-4 are derived using n = 20 and At=0.05 again for a single node in 
the middle of the tether. Notice no improvement on the error between the two results. Lastly, Figures 

2.2. 1- 5 and 2.2. 1-6 are shown using n=100 and zl/=0.05 for the middle node of the tether. Notice a 
remarkable improvement in the reduction of error between the exact string solution and the tether 
finite difference code. 
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Figure 2.2. 1-1: Comparsion of String Solution of Wave Equation and FD Code 

(20 nodes, 0.1 sec time step) 
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Figure 2.2. 1-2: Difference in Displacement of String Solution of Wave Equation and FD Code 

(1 node) 
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Figure 2.2. 1-3: Comparsion of String Solution of Wave Equation and FD Code 

(20 nodes, .05 sec time step) 
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Figure 2.2. 1-4: Difference in Displacement of String Solution of Wave Equation and FD Code 

(1 node) 
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Figure 2.2. 1-5: Comparsion of String Solution of Wave Equation and FD Code 

(100 nodes, .05 sec time step) 
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Figure 2.2. 1-6: Difference in Displacement of String Solution of Wave Equation and FD Code 

(1 node) 
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2.2.2 Case 2 


Case 2 consists of a series of convergence tests perfonned on a single algorithm to verify the stability 
of a given result. The convergence tests will also be used to indicate the step size in space and time 
that can be used for a given result. 

Results 

First, convergence on spatial step size was considered. Multiple simulations were performed with 
decreasing step size (increasing number of nodes) using the 3-point finite difference fonnulation and a 
fourth-order Runge Kutta integration routine. Overall tether length was considered as a metric for 
comparison, determined as the sum of the strained tether segments. Three cases are compared, n = 50, 
100 and 150 with the criteria, c=At/Ax >=1 ensured for each case. These results are shown in the 
following 2 Figures. Figures 2. 2.2-1 and 2.2.2-la shows plots all three cases over one orbit and at the 
end of the orbit. Notice that converged results are observed in the 100 node case. Figures 2.2. 1-2 and 
2.2.1-2a plot the difference in the 50 and 100 node from the 150 node case demonstrating sub .15 m 
tolerance in convergence for the 100 node case. 



Figure 2.2.2-1: Nodal Convergence for 3 Point FD Formulation 
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tether length (avg, m) 


Nodal convergence for 3 pt. FD formulation 



Figure 2.2.2-la: Nodal Convergence for 3 Point FD Formulation 



Figure 2.2.2-2: 


Difference in Nodal Convergence for 3 Point FD Formulation 
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Figure 2.2.2-2a: Difference in Nodal Convergence for 3 Point FD Formulation 


Summary 

Convergence of both the spatial step size and time step size has been successfully performed. The 
results demonstrate that for a 90 km tether, 100 nodes and . 1 second time step provides desirable and 
accurate performance for this formulation. 

2.2.3 Case 3 

Case 3 consists of validation between two different integration routines in the tether finite difference 
model. The case to be validated against is the standard tether model consisting of 3 point finite 
difference and 4 th order Runge-Kutta integration routines. The case validated is one using the same 3 
point finite difference routine but an Adams-Bashforth-Moulton 4 th order integration routine is 
employed. 

Results 

Both models were compared using n = 20 and At=QA and recorded for a complete orbit. Figures 2.2.3- 1 
and 2. 2. 3-2 show the results of the comparison between the two models when observed for a single 
node (middle node). Figures 2. 2. 3-2 is the error difference between the two models for the middle 
node over the entire orbital period. 
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Figure 2. 2.3-1: Comparison of 4 th Order RK and ABM FD Code 


(1 node) 


Difference in displacement of 4th order RK and 4th order ABM FD code (one node) 



Figure 2.2.3-2: Difference in Displacement of 4th Order RK and ABM FD Code 

(1 node) 


2.2.4 Case 4 


Case 4 consists of comparison between two different spatial differentiation routines in the tether finite 
difference model. The case to be validated against is the standard tether model consisting of 3 point 
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finite difference and 4 th order Runge-Kutta integration routines. The comparison case employs a 5 
point finite difference routine but a 4 th order Runge-Kutta integration routine. 

Results 

Both models were compared using n = 20 and At=0 . 1 and recorded for a complete orbit. Figures 2.2.4- 1 
and 2. 2.4-2 show the results of the comparison between the two models when observed for a single 
node (middle node). Figures 2. 2. 4-2 is the error difference between the two models for the middle 
node over the entire orbital period. 


Comparison of 3 point and 5 point FD code (one node) 



Figure 2.2.4-1: Comparison of 3 and 5 Point FD Code 

(1 node) 

Difference in displacement of 3 point and 5 point FD code (one node) 
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Figure 2.2.4-2: Difference in Displacement of 3 and 5 Point FD Code 

(1 node) 
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2.3 Validation Results for Type-B Cases 

The type-B validation cases compare results across multiple code platforms and looks for 
commonalities and differences. While not a definitive means of validation, corroboration of important 
dynamic phenomena by several codes is desired from the standpoint of code validation. Note that this 
cross-comparison may be used in the future to validate more efficient coding procedures, such as those 
using an assumed-modes method for example. A summary of this work has been prepared in formal 
paper format and is included at the end of this section. 

The following three primary cases are first considered. These cases are chosen to isolate potential 
dynamic effects. Each of the following cases consider a uniform tether without bending, tether mass of 
.0048 kg/m. 

Data for comparison: 

1) total tether length 

2) transverse vibration at nodes 6 and 12 (quarter and mid node) 

Case 5a) 

Brief description: massless, symmetric tether rotating in free space 
Purpose: Observe effects of rotation on a symmetric tether (no end masses) 

Parameters: mi= 0kg, m 2 = 0kg, 1 = 10,000 m, 

Initial unifonn stretch of .0002, rotation rate = .01 rad/s 
Number of nodes =23 

No damping, bending, aerodynamic, gravity or other effects considered 
Run time: up to 5000 s 

Response for this case is shown in average tether length (Figure 2.3-1) and transverse vibration at 
mid and quarter nodes (Figure 2.3-2) 

Case 5b) 

Brief description: Massive, symmetric tether rotating in free space 
Purpose: Observe effects of rotation on a symmetric tether 

Parameters: mp 1000kg, m 2 = 1000kg, 1 = 10,000 m, 

Initial uniform stretch of .0002, rotation rate = .01 rad/s 
Number of nodes =23 

No damping, bending, aerodynamic, gravity or other effects considered 
Run time: up to 5000 s 

Response for this case is shown in average tether length (Figure 2.3-3) and transverse vibration at 
mid and quarter nodes (Figure 2.3-4) 

Case 5c) 

Brief description: Non-symmetric tether rotating in free space 
Purpose: Observe effects of rotation on a non-symmetric tether 

Parameters: mp 1000kg, m 2 = 500kg, 1 = 10,000 m, 

Initial uniform stretch of .0004, rotation rate = .01 rad/s 
Number of nodes =23 
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No damping, bending, aerodynamic, gravity or other effects considered 
Run time: up to 5000 s 

Response for this case is shown in average tether length (Figure 2.3-5) and transverse vibration at 
mid and quarter nodes (Figure 2.3-6) 

Case 6) 

Brief description: Non-symmetric, fully-rotating tether in earth orbit 
Purpose: Observe tether behavior with gravitational effects 

Parameters: mi= 1000kg, m 2 = 500kg, 1 = 10,000 m, 

Initial uniform stretch, stretch ratio = .0004, rotation rate = .01 rad/s 
Number of nodes =23 

No damping, bending, aerodynamic or other effects considered 

Earth orbit: r p = 6723 km, r a = 16373km, e = .4178, spherical-earth gravity model 

Run time: up to 1 orbit (12.6 ks) 

Response for this case is shown in the motion of the center of mass (Figure 2.3-5), average tether 
length (Figure 2.3-6) and transverse vibration at mid and quarter nodes (Figure 2.3-7) 

TTU simulation results from case 5a: 
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Figure 2.3-1: Average Tether Length over Case 5a 
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case 1 , Transverse vibrtaion at quarter and mid nodes 
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Figure 2.3-2: Mid and Quarter Node Vibration, Case 5a 


TTU simulation results from case 5b: 


Aug. length of tether: m1=m2=1e3kg, I=1e4km, n=23 
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Figure 2.3-3: Average Tether Length over Case 5b 
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Transverse vib. of tether node 6, 12: m1=m2=1e3kg, I=1e4km, n=23 



TTU simulation results from case 5c: 


x -| q 4 Awg. length of tether: m1=1000kg, m2=500, l=10km, n=23 
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Transverse vib. of tether node 6, 12: m1=1e3kg m2=.5e3kg, 1=10 km, n=23 



TTU simulation results from case 6: 

non-symmetric tether, ml= 1000kg, m2= 500kg, l=10km, initial stretch=.0004 
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Figure 2.3-7: Orbit of Tether Center of Mass, Case 6 


x -|0 7 orbit of CM - symmetric 1000kg, 10km tether 
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x -|q 4 average length of tether - 1000kg/500kg, 10km, earth orbit 



Trans\erse vibration of mid and quarter node - 1000kq/500kq, 10km, earth orbit 
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Figure 2.3-9: Mid and Quarter Node Vibration, Case 6 
2.4 Discussion of Results 

Sections 2.2 through 2.3 provide a variety of comparison of simulated tether response with kn own 
response functions as well as side-by-side comparison of tether response as predicted by independently 
developed models. These models are based on an identical set of equations of motion, but vary matters 
such as base coordinates and propagation approaches. Reasonable measures were taken to help 
minimize numerical error, such as the use of adaptive time step integration routines and monitoring the 
system energy (no damping is considered in these models). The examples are arranged in such a way 
that some cases consider tether motion absent of external excitation effects, while others consider the 
tether in a gravitational field (spherical gravity model). In general, these comparisons taken together 
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provide a good level of verification for the developed model, leading to a validated model for use in 
simulating the MXER tether system. Some observations are made on some of the specific cases 
considered here. In particular comparisons between a finite-difference based model (described in 
section 1) and an assumed-modes model (developed at the University of Delaware) will be discussed. 

Case 5a considers a massless tether with symmetric endmasses, rotating in a deep-space environment 
(no external effects). Figures 2.3-1 and 2 show the predicted response for this case. Since the tether 
is massless, it is expected to behave as a linear spring vibrating in an extensive fashion according to 
the initial conditions. The models correspond readily and predict primary longitudinal period of 
approximately 13 seconds. 

Case 5b extends case 5a by allowing dynamic effects of the tether mass (massive tether considered). 
Figures 2.3-3 and 4 show the predicted response for this case and again the models in general show 
positive agreement. Again, amplitude and period of longitudinal vibrations correspond (8 meters peak 
to peak with 15 second period). However, a departure in the predicted results is also observed, the 
finite-difference model developed in section 1 predicts a beating phenomena to occur that suggests a 
coupling between the longitudinal and transverse modes of vibration. Such a coupling is not seen in 
the assumed-modes-based model. 

Case 5c considers a massive tether with non-symmetric end-masses, again rotating in a deep space 
environment, with results shown in Figures 2.3-5 and 6. As in case 5b, the models correspond in some 
aspects, such as amplitude of longitudinal vibration and period, but depart again in the suggested 
beating phenomena. The finite difference approach seems to predict a coupling in energy between the 
transverse and longitudinal vibration modes. 

Finally, case 6 adds external forces to the model of case 5c, by placing this massive, non-symmetric 
end-mass tether in a circular earth orbit. The primary effects of longitudinal and transverse vibration 
are shown in Figures 2.3-7 and 8 for the two models. These primary effects are consistent in terms of 
amplitude and frequency data. Coupling between these modes is not evident in either approach. 

In general, the comparison between the two alternative spatial integration approaches demonstrate 
reasonable agreement in most expected, linear-type behavior. The agreement in response is near 
complete in cases 5a and 6 considered. In the two central cases (cases 5b, c), considering a massive 
tether without external excitation, a phenomena (herein tenned beating) is observed in the finite- 
difference-based approach, but not in the assumed-modes method. This dynamic behavior seems to be 
described as a coupling between the transverse and longitudinal modes, such that vibrational energy is 
passed from one to the other. This takes place with no net change in energy over time, as expected in 
these cases that do not consider any type of damping. The explanation, source or even verification that 
this is a dynamic phenomena and not numerical is currently unavailable. On one hand, it is possible 
that this is a nonlinear effect that is seen in the finite difference approximating technique, but not 
spanned by the modes used in this analysis. Alternatively, the approximations incurred by the finite 
difference method could generate unexpected numerical behavior. Due to correspondence in case 5 a, 
the former is currently expected and an effort to identify the contributing nonlinear dynamic effect is 
underway. It should be noted that this coupling behavior is not observed in case four for an earth- 
orbiting tether. It is proposed that the external effects of gravity shift the frequencies in such a way 
that the beating phenomena is not seen in the time span considered, or may be eliminated completely. 

It should also be noted that in case 6, agreement is demonstrated in the tether modal behavior. 
However, it is noted that observing the motion of the center of mass would predict two different 
phenomena. In the finite difference approach, the center of mass is seen to migrate above and below 
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the keplerian orbit, as expected due to the gravity gradient effects on the rotating massive tether. 
However, due to the polynomial modes selected, the assumed-modes method does not permit motion 
of the tether center of mass off the local reference coordinate system, imposing an arbitrary constraint 
on its overall motion. 

3.0 BASIC IMPLEMENTATION AND DEMONSTRATION 

(In Selected Software Platform) 

3.1 Tether Dynamic Models 

(F ormulation- Specific Algorithms) 

The dynamic models and specific formulations appropriate for the MXER tether system have been 
demonstrated in Section B 1.0. The results of that section demonstrated a number of formulations, all 
of which require computational treatment (most through numerical solution techniques). This section 
(Section B, 3.0) will demonstrate specific algorithmic routines in diagram form through the use of flow 
charts. The algorithm flowcharts will first be presented in the most general form, showing the primary 
functions that need to be implemented and their relation to program flow control. Subsequent 
flowcharts will demonstrate the specific processes performed in the each function, and will present the 
variety of functional approaches dictated by the various formulations. A few notes are made for the 
following algorithm flow charts. The algorithms as suggested here will be implemented with strong 
dependence on the use of functions. Information passed to and from these functions will be presented 
within the flowchart. Further, unless additionally specified, the scope of any defined variables will be 
limited to its own functional block. Functional blocks will be demonstrated in the following 
flowcharts with the use of dashed border lines. An additional note regarding variables is provided. To 
add detail to the function-specific flowcharts, some variable names will be introduced and used in a 
general fashion. These variable names are introduced as a guideline rather than a necessity for 
developing subsequent code. Of more importance is the variable structure suggested in these 
flowcharts, with the tether state contained as a one dimensional array, S containing specific state 
vectors in the order of r, v, %. As a final note, this section is designed to detail the computational 
implementation and ease code design. To do this, the flow charts will be grouped together at the end 
of the descriptive text, with each flow chart occupying one page allowing the flowchart to be printed or 
copied and tagged as a description to each computational function. 

1 . The general program formulation for evaluation tether dynamics over a period time is first presented 
as shown in Figure 3.1-1. 

Functions to define initial conditions as well as techniques for post-processing are general to most 
problem formulations and will be presented later in this document. 

2. The first specific instance of the general formulation of the dynamics main code is presented next. 
This first formulation implements a fourth-order Runge Kutta time integration routine and a three- 
point finite difference spatial approximation routine as shown in Figure 3.1-2. Function-specific 
details then follow in flowcharts Figure 3.1-3 for the time integrator and Figure 3.1-4 for the 
EOM/finite difference routine. 

Note that the evaluation of forces on the tether still exists as a function block. Evaluation of the tether 
force merits significant discussion and detail which will follow later. Note that full consideration of 
tether forces will most likely include a number of dynamic, gravitation and aerodynamic models which 
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will be called upon as independent functions made available from a variety of sources (indicated in a 
general format in Figure 3.1-5). Current implementation however considers only a simple gravity 
model that is embedded completely within the EOM function. 

3. A second implementation of the general tether dynamic simulation routine considered will 
implement an alternative integration routine based on the Adams-Bashforth-Moulton technique 
(Hoffman, 2000). This implementation will be based on the same Dynamics-Main function 
demonstrated in Figure 3.1-2 with the new (ABM) integration technique inserted. The detailed 
implementation of the ABM integration routine is demonstrated in Figure 3.1-6 

The general algorithm determining initial conditions is given in Figure 3.1-7. This code takes in tether 
design and orbit parameters and returns the tether initial state as well as the properties structure. 



1 

( M ) 


Figure 3.1-1: Flow Chart, General Tether Dynamics Integration 
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Figure 3.1-2: Flow Chart, Overview of Runge Kutta Routine 
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Figure 3.1-3: Flow Chart, Details of 4th Order Runge Kutta 
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EOM/Finite difference routine 
3-point 


Extract r, v, % 
rx=S(l :3:3*n) 


Evaluate % - Central difference: 
for i=2:n-l 

X'(0 = (x( i+1 )- xO" 1 ))/ 2 ^ 

for i=l 

X'(0 = (x(' +1 )- X(0)/Ax 
for i=n 

X'(0 = (x(0- x('-l)V Ax 

v y 


Evaluate r|, t|' 

for i=l :n 

Hi - (Xxi + Xyi + Xzi ) 
H i - (x xi + X yi + X zi ) 


Evaluate F 

for i=2:n-l 

F =f[p.S) 

(force per length) 

for i=l,n 

F=/(/77,S) 

(force) 


Evaluate v' - Central difference: 
for i=2:n-l 

v'(i) = (v'(i + l)- v'(i-l))/2Ax 
for i=l 

v'(i) = (v'(i+l)- v'(i))/Ax 
for i=n 

v'(i) = (v'(i)- v'(i-l))/Ax 

v y 


Solve EOM's: 

for i=2:n-l 

Av k = At/p * EA(ri(r|-l)x'k+ if Xk) n' 2 + F k 
for i=l,n 

Av k (i) = At/m* EA((ti-1)xk+ F k 
for i=l :n 

Ar = Atv 
Ax - At v' 


Rebuild dS/dt 
AS(l:3:3*n)= Ar x 

I dS/dt 


Figure 3.1-4: Flow Chart, Overview of 3-Point Spatial Difference Routine 












Figure 3.1-5: Flow Chart, Obtaining Tether Force Components 



Figure 3.1-6: Flow Chart, Adams-Bashfort Mouton Integration 
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Initial Condition Routine 



Figure 3.1-7: Flow Chart, Obtain Initial Conditions 


4.0 BASIC APPLICATION OF THE DYNAMIC MODEL 

(As A Simulation, Analysis and Design Tool) 

Based on the models (section 1) and algorithmic implementation (section 3), general software is 
developed for implementation. This implementation has been conducted on three primary platforms, 
Matlab, C++ and Adams. Code listings for the developed software are included in Appendix B.4. 
This code has been demonstrated as an analysis and design tool in section A. 1 of this report for 
example, where the dynamic modeling code is used as a basis for stochastic sampling of the tether 
propagation process. 
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5.0 DEVELOP INITIAL AND BOUNDARY CONDITION MODELS 


Initial conditions for the tether are derived from a tether model that is reduced to the level of analytical 
solution. The assumptions that give this reduction include a uniform tether, constant rotation rate and 
absence of gravity gradients. The initial conditions are derived here with a reference location selected 
as the instant center of the tether at the initial state as shown in Figure 5.0-1. 



Figure 5.0-1: Schematic of Tether Configuration to Provide Steady State Intial Conditions 


Consider nominal tether motion with instant center of rotation as shown in Figure 5.0-1 . The initial 
parameters given are; end masses, tether nominal length, unifonn tether, end mass and center of mass 
velocities. Additional parameters in Figure 5.0-1 are determined as; 

Y = 0 + (/) 6-1 


and 



6-2a-c 


6-3 


is the acceleration of the instant center in a direction aligned with the tether. The tether’s initial state 
referenced from the instant center is shown in Figure 5.0-2. 
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dr 
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pAp/r 2 dr * 
pAdr(y 2 r+a 0 ) 

Figure 5.0-2: Free Body Diagram of Tether Element 
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Equations of motion for an element in the initial tether is; 

N ' - P A = ~P A (r 2 r + a 0 ) 


with 


Equation 6-4 becomes 


with the solution, 


N = EAs = EA(r' - l) 
N' = EAr" 


r- + PE r =P' 


7 


P/ 2 -a, 

\/ r p J 


= C x cos(/k) + C 2 sin(/?x) + 


f / A 
V/ 7 


>9 = 


/w 2 


6-4 

6-5a,b 

6-6 

6-7 

6-8 


The two constants in Eq. 6-7 are detennined from the following two boundary conditions. 
At x=a, 

m, ' ' A 


''(a) = —P 
V EA 


r 

v/ /> 


-(;f 2 tf + tf 0 ) +1 = -Cyj3sm(j3a)+C 2 j3cos(j3a) 6-9 


At x=h, 


r '^ =, ^EA ^ 2b + a °} 


- P/ 


J 

\ 

' p J 


+ 1 = -Cyj3 sin {fib) + C 2 p cos (fib) 


6-10 


Yielding the coefficient solution. 
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+ 1 
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'p J 


6-11 


With coefficients Ci and C? known, the initial state is given in Eq. 6-7 for the position coordinates and 

r' = -C x Psm(f3x)+ C 2 /3cos(fix). 6-12 


Figure 5.0-3 below plots a sample resulting initial value ( 77 ) for a tether case. 
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Figure 5.0-3: Tether Strain Over Length for Steady State Initial Condition 
6.0 RESULTS AND CONCLUSIONS: 

This report has described the significant activities that were performed by the TTU organization during 
2004-2005 on the MXER tether system. Significant results of these activities include the following; 

1 . A capture mechanism was developed, designed, fabricated and tested. This mechanism 
demonstrated robust and reliable capture within the design parameters of the tests. Further, 
within the assumptions of tether motion prediction capabilities, it is believed that this 
mechanism will provide allow the MXER system to achieve its necessary rendezvous and 
capture process. While the final stage of capture, release, was not specifically addressed, it is 
considered to be a tractable issue within the design capabilities of the proposed mechanism. 

2. A series of other mechanisms were proposed and are considered realistic alternatives to 
achieving MXER tether system capture. The advantage of multiple options is believed to be 
significant. Further development, design and testing would be necessary to achieve the same 
confidence in these mechanisms as in the primary proposed mechanism (Quad-Trap). 

3. A facility was developed that provides a dynamically realistic environment for testing scaled 
MXER tether capture mechanisms. 

4. Dynamic models suitable for the elastic, rotating MXER tether have been developed, along 
with associated algorithms. Multiple fonnulations of these models are demonstrated. 

5. A rigorous validation of the dynamic modeling capabilities for the MXER tether system was 
conducted. This process consisted of several components, including comparison with closed 
form solutions, limiting cases and comparisons over multiple codes developed by different 
groups. The result includes validation of specific dynamic models as well as the development 
of a series of test cases. 
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6. Significant nonlinear effects are recognized in the tether behavior that have not been reported 
in the literature. The ability to accurately characterize such nonlinear behavior represents the 
level of fidelity that may be desired in such code. 

7. The models and algorithms developed in this work are recognized to be limited in tenns of 
processor demands relative to other algorithms that have been developed. However, the 
primary, significant contribution of this effort is to serve as a tool that will allow the 
development and verification of time-efficient codes. 
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SECTION C— ERROR SPACE PREDICTION 


This section presents a formal paper prepared for submission to AIAA journal transactions 

Error Space Prediction and Dynamic Response of 
Passive Capture Mechanisms in Tether Momentum Exchange 


Momentum-exchange/electrodynamic reboost (MXER) tether systems may enable high-energy missions to the Moon, 
Mars and beyond by serving as an "upper stage in space". Existing rockets that use a MXER tether station could 
double their capability to launch communications satellites and help improve US competitiveness. A MXER tether 
station would boost spacecraft from low Earth orbit to a high-energy orbit quickly, like a high-thrust rocket. Then, 
using the same principles that make an electric motor work, it would slowly rebuild its orbital momentum by pushing 
against the Earth's magnetic lield-without using any propellant. One of the significant challenges in developing a 
momentum-exchange/electrodynamic reboost tether system is in the analysis and design of the capture mechanism 
and its effects on the overall dynamics of the system (Sorensen, 2001) 1 . This paper will provide qualitative and 
quantitative information about the nature of the capture window due to parametric uncertainties in the MXER 
tether and will consider the dynamic effects of capture at various locations within this error window. 


Nomenclature 
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A, I 
E 

nil 

m 2 

v 

r 


e, 77 

F 

X 

K 
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R 

P 

C h C 2 , p 
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Pxt Iyy> fvv 


CT 


tether material density 

tether cross sectional area, area moment of inertia 

tether modulus 

facility mass 

capture-end mass 

velocity 

position 

strain 

force 

spatial derivative of position 
absolute value of % 
state vector 
acceleration 
orbital radius 
gravitational constant 
constants of integration 
centroid 

mass moment of inertia 
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Introduction 


A critical component in the development of a MXER tether system is an appropriate capture device, 
which will accommodate kinematic and dynamic variation between the tether and payload during the 
rendezvous and capture process. Further, it is desirable to do this as efficiently as possible through the 
design of a passive capture mechanism. Therefore, a necessary step in the design of this capture 
mechanism is to define the kinematic and dynamic requirements for this device. This paper will attempt 
to address this need through developing qualitative and quantitative infonnation about the nature of the 
capture window due to uncertainty in input parameters in the MXER tether model. This data will fonn a 
capture window and will be used to guide the design of an appropriate capture mechanism devices to 
ensure reliable rendezvous between the MXER tether system and its target payload. This information 
will consist of an estimation of the geometric nature of the capture window, both in shape, size and 
orientation relative to the tether as seen by the payload. The capture window will be based on a heuristic 
assumption for its geometric form. Confinnation and quantitative evaluation of the properties of this 
error geometry will be determined stochastically through the use of a Monte Carlo simulation based on 
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uncertainty in tether parameters within an elastic numerical tether model. Transformations are 
introduced to express the stochastic error data in terms of geometric parameters relative to the tether 
motion. Infonnation from these analyses can be used in both capture mechanism design and will indicate 
bounds for the capture process based on stability of tether motion. 

Past research in the field of space tether modeling has seen a significant amount of work in modeling the 
dynamic characteristics of these systems (see for example [1-3] (Misra, A.K. and Modi, Beletsky and 
Levin, Quadrelli, B.M., and Lorenzini, E.C)). The proposed space tether application as a momentum 
exchange vehicle requires infonnation about the location of the tether tip at time of capture and the 
sensitivity of the tether tip location to parametric variance. However, studies of sensitivity of tether 
model dynamics due to parametric variance are much more limited in the field of space tethers. The few 
examples that are available are briefly mentioned. Tragessar et ah, consider stochastic error effects on 
tether-based aerocapture [4]. Pelaez and Lorenzini consider sensitivity of unmotorized deployers such as 
the SEDS (Small Expendable Deployment System) on tether-based satellite injection [5]. Finally, Leamy 
et al. develop sensitivity coefficients for electrodynamic response of bare- wire tether systems [6]. While 
these papers develop considerations based on tether sensitivities, a sensitivity model that will allow 
extension to the rendezvous capture process associated with the MXER tether system is needed and is the 
focus of this paper. 

The remainder of the paper will proceed in the following manner. First, a summary of the numerical 
tether model used for simulation is provided. Then, the Monte Carlos simulations are performed to result 
in predictions of the capture error box. Appropriate transformations are defined to express the error 
values into a geometric form for the capture window and are evaluated for the various simulations 
performed. Analysis of the dynamic behavior over these capture windows is conducted and presented 
through a series of phase portraits. Finally, the results of these efforts and their implications on capture 
mechanism design are presented. 


Overview of the Tether Momentum Exchange Model 


The equations of motion for an elastic space tether within its orbital plane are derived using Newton 
mechanics, in a manner similar to that presented by Kohler et al, (1978)". The resulting equations, which 
consider both axial and flexural stiffness of the tether, are given as, 


pAvj = Ea[ii'Xj + {ij - \)ijx'j 1 + E1[tj{k"x k + k'x[ ) - 2 ij'k'x, \] 3 + F } , (1) 


(j,k) = (x,y), (y,x), 

with tether parameters p, E, A, and / representing tether density, modulus, area and area moment of 
inertia respectively, Fj the external force on the tether, trailing primes represent spatial derivatives (with 
respect to the geometric coordinate), overdots represent time derivatives and parameters r|,x and k, are 
defined as, 


r/ = s + 1 = 


ds * 
ds 


(2) 


and 


2T = 


dr 

ds 



K ~ yjx'x+z'y+z'l 

K'=\x:x:+z'x+x:z:y' 

+ x'x + lzu:+ zx +xk 2 


(3) 
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(4) 

(5) 

(6) 



with s the tether strain and s and s* the tether parameteric coordinate based on stretched and unstretched 
length respectively. The equations defined in terms of the time and spatial derivatives of r defined in the 
state vector, X as, 



The equations of motion (1) are second order in time and fourth order in space. For this system, four 
boundary conditions and the initial state are required. The boundary conditions will be defined as the 
end force and moments that satisfy the equations of motion for the end bodies, given in the following 
equations as, 


pAl\ 

= 

EA(rj- l)x 

(8) 

77 (o,0 

77 (ofi 


P A i 2 X 

m b2 a b2 

EA(ti- 1)jc 

F b2 

(9) 
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77 (/,) 


(o' 

II 

0 

(10) 

It' M = 0 

(11) 


with mti, in/, 2 , F/,/,and F 7,2 the mass and force at the tether boundaries. Initial conditions for the tether are 
derived by considering the tether in an instantaneous sense, with the tether lying along the axis of true 
anomaly and its motion described as pure rotation about an accelerating instant center. In this initial 
case, gravity gradient effects are also ignored. This allows the tether to be described with a partial 
differential equation that can be solved in closed form; 


EAr" + py r = p\ 



( 12 ) 


with r describing the initial position of points along the tether relative to the instant center, y 1 the 
rotation rate of the tether about its instant center, // the gravitational constant, R the distance of the tether 
center of mass (from earth center) and cio the acceleration of the center of mass. From Eq. (12), the 
initial state can be solved as, 



with x the undeformed positions of the tether a coefficients Cj and C? determined from the equations of 
motion for the tether end-bodies. 

The tether equations of motion given in Eq. (1) along with the tether along with boundary conditions 
initial conditions are solved in time to give a simulation of the tether. The numerical algorithm employed 
in this work uses a finite difference approach to transfonn the governing equations for the tether into a 
system of ordinary differential equations that are integrated in time using a This model used a Runge- 
Kutta-Fehlberg formula [7]. The solution is validated through comparison with solutions presented by 
Fried [8] and Henghold and Russell [9] for an earth-based dynamic cable system, and through 
comparison with exact derivations of a space tether modeled as a two-mass orbital system connected by 
an extensible mass-less rod. 
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Monte Carlo Simulation Approach 

A prediction of the capture window as seen by an observer on the payload is constructed based on 
uncertainties in the tether parameters. The capture window will be defined as the spatial region of 
computed tether tip positions relative to the payload at the point of capture. For this effort, the capture 
window will be initially described geometrically as an ellipsoid, with an assumed normal distribution 
(i.e., the relative position or error between tether tip and payload is distributed normally). This window 
will be called an error ellipsoid. The spatial region of tether tip error will be generated using the Monte 
Carlo technique and will consider stochastic uncertainties in the tether input parameters. In particular, 
the Monte Carlo method employed here is that which Morgenthaler [10] has tenned model sampling with 
the primary variability in the system assumed to come from uncertainty in tether parameters. 

The Monte Carlo simulation will be conducted as follows. Numerical simulation of the MXER tether 
model is conducted using the model given in part II, starting from a steady state condition of the tether on 
orbit and ending when the tether reaches the point of capture. This numerical simulation is conducted 
repetitively with stochastic variation in tether parameters in place to yield variability between the payload 
and tether tip at the defined point of capture. The collection of this output positional variability then 
results in the data for the desired capture window. With an eye toward the eventual purpose of the error 
window analysis, the capture window will be assumed to take the geometric fonn of an ellipse (ellipsoid) 
which can be described in terms of center and major and minor axes. The span of the capture window 
can now be defined based on the major and minor axes of the error ellipse. The distribution of this data 
within the ellipse will also be assumed to be normally distributed. A comparison of the Monte Carlo 
simulation results with this assumed form will be conducted to ensure the reasonable nature of these 
assumptions. 

The Monte Carlo method is selected for use in this paper based on its generality and adaptability in 
implementation. As indicated by Kushner et ah, [11] the Monte Carlo method is advantageous in highly 
nonlinear problems, is insensitive to system dimensionality and is demonstrated to provide a robust result 
for sensitivity estimates. Further, the programming algorithm based on the Monte Carlo approach 
provides a means to avoid numerical integration or derivation and provides reasonable convergence as 
demonstrated by Thompson and Cluett [12], [13]. 

Implementation 

The primary details for implementing the Monte Carlo simulation involve definition of parameter 
variability and proposed models for defining the point of capture. Stochastic variability in the tether 
parameters is assumed to be purely random or white noise with no correlation between parameters or 
successive parameter values [13] centered about the parameter nominal value with a normal distribution. 
Table 1 lists the specific tether parameters that will be varied along with their nominal (mean) value and 
standard deviation. The mean parameter value is assigned as the baseline nominal parameter value while 
the standard deviation is assigned and 1/3 the maximum variation expected of a given parameter. Table 2 
provides the remaining parameters used in defining the tether and payload system model. 
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Table 1: Specific Tether Parameters 


Parameter 

% Variation 

Capture Mechanism Mass 

0.01 

Ballast Mass 

0.01 

Total Tether Length 

0.01 

Cross Sectional Area 

1 

Elastic Modulus 

2 

Mass Moment of Inertia 

0.01 

Mass Per Unit Length 

0.01 


Table 2: Tether and Payload Parameters 


Constants 


1.1. 1.1. 1.1.1 Tether 

Facility 
Orbit 



3.986E+14 

Perigee (m) 

6,785,136.00 

Payload Mass (kg) 

1000, 2500 

Apogee (m) 

14,823,136.00 

Capture Mass (kg) 

250 

Time for one orbit (s) 

11,176.25 

Ballast Mass (kg) 

11000 

SMA 

10,804,136.00 

Tether Length (m) 

90000 

e 

0.371987172 

Cross Sect Area (rrf ) 

1.9635E-05 

Perigee Velocity (m/s) 

8977.690593 

P (kg/m 3 ) 

970 

Apogee Velocity (m/s) 

4109.444293 

Payload orbit 


co of tether fac. (rad/s) 

-0.013881192 

Orbit Radius (m) 

6694275.739 



Vp (m/s) 

7716.44186 



Time for one orbit (s) 

5450.876936 








The overall process consists of repetitive evaluations of individual tether simulation. At the beginning of 
each tether simulation, values for in the input parameters with variance are selected in a random fashion. 
Each tether simulation will begin with the tether in a state of unifonn strain (0.4% strain) one orbit prior 
to capture. The initial state will assume to be a known condition and therefore will not be considered as a 
parameter of uncertainty in this analysis. The tether simulation will continue to the point of capture, 
which is defined below. Upon reaching the point of capture, the error in capture (position of tether 
relative to payload) is determined and stored and the individual tether simulation is complete. The 
procedure will then continue with random variation in input parameters defined. The stopping criterion 
is given as a pre-defined number of iterations. 

Defining point of capture: 

The resultant error in capture that occurs in each tether simulation is recorded as the vector from the 
payload to the tether tip at the point of capture. Since the goal of this work is to develop a capture 
window as the starting point for the design of a passive capture mechanism, point of capture will be 
considered the ideal capture position that occurs (i.e., point of minimum distance between tether tip and 
payload). Note that the ideal capture position is a function of the orbits of the tether and payload. 
However, the scope of this paper considers variability in the tether model only and therefore will 
consider the payload orbit in an ideal sense. It still remains that the capture process is a function of space 
and time and could be specified in varying ways. This paper will consider the capture event over three 
possible definitions for point of capture; the first temporal, the second spatial, and the third a 
combination of space and time. 
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Capture method 1 : In the first method, the point of capture is defined to occur at a fixed, specific point in 
time. Capture error is measured as the position vector from the payload to the tether tip at the point of 
capture. The location of the payload is fixed in space and defined such that zero capture error results 
when the nominal tether parameters are used. Since the simulation will nominally consider one full orbit 
of the MXER tether, the end of period of the tether will be chosen as the time for capture, and this 
method will be referred to as End Of Period (EOP). 

Capture method 2: In the second method, the point of capture is defined to occur when the tether tip 
reaches its minimum distance to the payload where the payload position is fixed and defined uniquely as 
the location of the tether tip at a specific point in time such that zero capture error results when the 
nominal tether parameters are used.. Since this method looks for the closest approach of the tether to a 
fixed payload point, it will be referred to as Closest Approach (CAP). 

Capture method 3: In the third method, the point of capture is defined to occur when the tether tip 
reaches its minimum distance to the payload where the payload is now defined as a point mass traveling 
on a keplerian orbit such that zero capture error results when the nominal tether parameters are used. 
Since this method considers an orbiting payload, it will be referred to as ORBit (ORB). 


Capture Window Results 

The resulting data from the Monte Carlo simulations over the three methods of capture are shown in the 
following graphs. This data is plotted in a tether-tip frame with the x axis the line connecting the two 
ends of the tether, the y axis normal to x and lying in the orbital plane and the zero error position defining 
the origin of these axes. The resulting error data is shown in Fig. 1-3 for tether capture method 1-3 
respectively. The data in each case represents 250 data points. 
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Figure 1, case 1 


Error at Capture (ORB) 
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Figure 2, case 2 


Figure 3, case 3 


Figure 1-3: Tether tip uncertainty plots for case 1, 2 and 3 
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It is desired to mathematically describe these capture windows parametrically such that they can be 
defined in terms of a frame corresponding to axial and transverse displacements of the tether tip. A 
transformation is defined by considering error plots as a collection of error particles each with equivalent 
weighting. The density of error points resulting from the Monte Carlo simulation gives the fonn and 
density of this error body. This body can be described in terms of its centroid and equivalent moment of 
inertia, both given in base-frame coordinates, with centroid defined as; 

n n 

Y J x i dm l Y,yt dm i 

X C = M T c = — (14) 

C n s C yi V 7 

Yj dm i 

i = 1 

and body inertia tensor; 


Yj dm i 


i = 


2 », ~Y j x i y i dm i 

i = 1 i= 1 

n n 

Ys x iyi dm i Ys x i dm i 

i = 1 i = 1 

with dint representing a single error particle with position x„ v,. and having unit mass. 


(15) 


With the equivalent inertia tensor defined, the principle inertia and principle directions can be 
detennined. The principle direction will define the orientation of the capture window relative to the 
tether while the principle inertias will be used to define the direction of major influence of the capture 
window. 


The capture window will now be described in the form of an ellipse (called error ellipse). Assuming a 
uniform density (albeit unknown) of the ellipse , the principle inertias defined in terms of the major and 
minor axes (a, b ) are equated to the numerically derived values, resulting in closed form solutions for the 
ellipse axes as follows; 


16/ 


a = 


3 

max 




y = 16/ ” in 


;r 2 /.. 


(17) 


or as preferred, defined as ellipse eccentricity and radius at perigee. 

The error ellipse can now be described in tenns of its eccentricity, axes, a and b and its orientation 
relative to the tether, 6. This ellipse that contains 99.7% of the error points will be called the 3a error 
ellipsoid, such that standard deviation for the major and minor axes is defined as; 


a b 
<7 “ - 3 ,CJh ~ 3 ’ 


(18) 


The results for the three cases considered above are shown in Table 3 and in Fig. 2. The table reports the 
results for the equivalent principle inertia’s, orientation, and error ellipse major and minor axes for the 
three cases considered in tether parameter error. The figure demonstrates an error ellipsoid plotted over 
the error estimate points for case 1 . 


Table 3: Error Ell 


lipsoid Parameters for cases 1, 2 and 3 



Centroid 

(m) 

Imax (m ) 

Imin (m 4 ) 

0 (deg.) 

a(m) 

b (m) 

a a 

(m) 

CTb 

(m) 

CASE 

1 

(10.8, -.11) 

1.3007e6 

1.0941e3 

8.6e-3 

86.9 

2.52 

28.7 

0.83 

Case 2 

(37.1,0) 

5.2739e7 

1.4194e6 

-2.1e-3 

142.2 

23.3 

46.9 

7.69 
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Analysis of Tether Behavior through Rendezvous with Payload 

The resulting error windows observed for capture definitions of cases 1-3 are now considered. It is 
observed that the results from case 2 and 3 seem well suited to the error ellipse definition, each with 
major axis generally aligned with the tether. However, the error ellipse generated for case one does not 
seem to provide as good a descriptor of the error data, rather the error data seems to be aligned along an 
approximate sinusoidal curve. Further, in this case the major axis of this ellipse is generally orthogonal 
to the tether axis. To better explain the motion resulting from case 1, the end of period evaluation, a 
survey of tether tip motion is conducted for a variety of inputs, as shown in figure 4. This data for the 
error window for case 1 could be reconstructed from this tether tip data by finding the minimum length 
vector from the payload position (0,0) and the tether tip path. 

The second item that is considered from this motion is the general distribution of the data that 
constructed the error windows. It is believed that as the variance of the input parameters is increased, the 
size of the error ellipses will also increase, while retaining their same general fonn. This premise is 
evaluated through numerical analysis of the data, now considering three cases, the original case called the 
100% variance case, and two additional (one with 200 % variance and the third with 333% of the original 
variance. The results of these analyses are shown in the blue markers in Figs. 1-3. 

It can be seen that in each case, the error window is increased, and in general retains its original form. 

The one exception is that in case three, it is observed that as the input variance increases past a certain 
level, the error in the tether normal direction increases significantly. It is proposed that this is due to the 
fact that a significant chance in the tether arrival time occurs, and the initially dominant effect, tether 
axial vibration, is over-ridden by a second effect, tether rotation. 

Summary and Conclusions 

This paper has performed an analysis of the capture error window for a MXER tether system based on 
uncertainty in tether parameters within a numerical tether model. Due to the nature of tether behavior, a 
system consisting of axial and transverse vibrations, it is postulated that the error bounds resulting from 
uncertainty in tether behavior could be best described geometrically through an error ellipse or an error 
ellipsoid. Plots of the error distribution based on uncertainty in tether parameters are created using a 
Monte Carlo simulation of the tether numerical model for several scenarios of error prediction. The plots 
of these error distributions (Fig. 1 a-c) indeed depict an error profile that is spatially distributed along two 
primary axes much like the major and minor axis of the ellipse. In case 1 and 2 (closest approach to a 
fixed point in space), the distribution is generally normally distributed about the zero error position, with 
the density dropping off at the boundary of the errors. In case 3 (closest approach to an orbiting payload) 
the error distribution again falls within an ellipse type pattern, with the error points distributed largely 
along the boundary of the ellipse. The error of case 3 may be better described as a cycloid with the 
primary error resulting from timing issues between the payload and tether (these are not accounted for in 
this model). 

A transformation is defined to convert error data into an equivalent error ellipse as shown in Eqs. 14-18. 
The results for all three cases are shown in Table III. The cases are averaged to result in a current 
estimate for the error ellipse, with properties given in the last row of Table III. The major axis of this 
error ellipse is approximately 10 times greater than the minor axis, and is aligned with the axis of the 
tether. Thus, the major axis represents axial displacements of the tether while the minor axis represents 
transverse vibrations. 


132 



Tether response due to the capture process is observed through phase portraits of average tether axial 
strain and tether rotation. The tether response prior to capture and after capture is shown if Figs. 4-6. It 
can be noted that the capture process maintains the periodic nature of the tether response, but greatly 
increases the amplitudes of response in both longitudinal vibrations and in rotational rate. The response 
demonstrated in Figs. 4-6 is apparently equivalent, but represent capture at three distinct points within 
the error ellipse (minimum, median and maximum error). This result demonstrates that, over the error 
represented in this error ellipse, that the jerk forces introduced in capture greatly outweigh the potential 
impact forces in the region of error considered. 


133 



Phase Portrait on Tether rotation, at A 


Phase portrait on Tether Length, at A 
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Figure 4: Phase portraits at point A 



Phase portrait on Tether Length, at B 
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Figure 5: Phase portraits at point B 


Phase Portrait on Tether rotation, at B 



Phase portrait on Tether Length, at C 
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Figure 6: Phase portraits at point C 


Phase Portrait on Tether rotation, at C 



The results of this work will be used in the design of appropriate capture mechanisms for the MXER 
tether system. These results imply that an ideal capture mechanism would be able to account for error 
that is distributed about spatial axes with the error along the most significant axis (tether axis) 


134 


approximately ten times greater than errors in the other axes (transverse axes). By design the capture 
mechanism to the expected error can result in optimal configurations in terms of weight savings, 
reliability, and maximizing the capture window. Several areas of additional work are suggested in 
reaching improved estimates of the error window for capture. While this work only considers error 
due to uncertainties in tether parameters, many other sources of error may exist, for example in initial 
conditions or on-orbit properties. In addition, in the transfonnation process defined, it should be noted 
that the error distribution is generally normal while the inertias calculated for the ellipse are based on a 
uniform distribution. This could be corrected through an ellipse moment of inertia definition based on 
both ellipse geometry and mass distribution. 
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SECTION D— QUAD TRAP CAPTURE MECHANISM PROTOTYPE CAD 

MODEL 

This appendix presents the CAD model used to construct the prototype Quad trap capture mechanism. 
First, isometric views of the Quad Trap prototype cad model are provided demonstrating the 
mechanism in an open and closed position (Figures Appendix D-l, 2). Following that individual 
component drawings are provided for all components that are designed and manufactured specifically 
for the Quad trap prototype mechanism (Figures App.C-3 through Appendix D-10). The remaining 
materials used to construct the prototype are commercially available materials. Finally, a photograph 
of the fabricated Quad trap prototype is provided in Figure Appendix D-l 1. 



Figure 1: Isometric View of the Quad trap prototype in an open configuration 
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Figure 2: Isometric View of the Quad trap prototype in a closed configuration 
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Figure 3: Quad Trap Exploded view, One Corner 
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Figure 4: Quad Trap Exploded View 
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Figure 5: Component drawing: Solenoid Mount 
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Figure 6: Component drawing: Bearing Tube 
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Figure 7: Component drawing: Corner angle 
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Figure 8: Component drawing: Capture Plate 1 





Figure 9: Component drawing: Capture Plate 2 



Figure 10: Component drawing: End Cap 



Figure 11: Quad Trap Prototype as Fabricated 
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